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Abstract 



This dissertation addresses the question of how to detect light top squarks at the 
upgraded Fermilab Tevatron collider. After a brief introduction to supersymmetry, 
the basic phenomenology of the light stop is reviewed and the current experimen- 
tal situation is surveyed. The analysis presented here is based on collider event 
simulations. The main decay modes accessible to the Tevatron are studied, feasible 
discovery channels are identified, and recipes for experimental analysis are proposed. 
It is found that stops with masses up to the top quark mass are liable to detection 
under these schemes with the data from a few years' running at the upgraded Teva- 
tron. With such an extended run, significant portions of parameter space may be 
probed. 
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Chapter 1 

Introduction to Supersymmetry 



This thesis investigates strategies for detecting a top squark (also referred to as the 
scalar top or, for short, stop) at the Fermilab Tevatron pp collider. The top squark 
is a representative of a new class of particles in the theory of supersymmetry. It is 
thought by some that the stop could be one of the most accessible particles of this 
class, and experimental physicists have invested a great deal of effort in the search 
for this particle. 

We begin in this Chapter with an introduction to supersymmetry and a 
brief exposition of its attractions for particle physicists. We characterize the class of 
practical supersymmetry models and offer a synopsis of the phenomenology relevant 
to our study. Then we introduce the top squark and discuss its properties in Chapter 
H We present the basic formulae and indicate how the particle may be light, and 
also why it might be expected to be light. We summarize the production and 
decay modes of the stop, and include a precis of current experimental results on the 
stop and other supersymmetric particles important to our work. In Chapter ||, we 
describe calculational techniques and the methodology of our analysis. 
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Our analysis of how experiments at the Tevatron may search for the stop 
begins in Chapter || with an investigation of the decay mode t\ — ► cZ\. We review 
the event topology and identify physics background processes which can mask its 
signal. We look in several channels, and by examining key kinematical features of 
the stop events we are led to propose a program of experimental procedures designed 
to maximize the discovery potential of the stop decays via this mode. We perform 
a similar study for the decay mode t\ — > bW\ in Chapter [|. Our conclusions and 
some summary remarks on the prospects for uncovering the stop at the Tevatron 
are presented in Chapter ||. 

1.1 Supersymmetry as a fundamental symmetry of na- 
ture 

The body of data presently available from investigations into the fundamental con- 
stituents of matter is consistent with a theoretical framework known as the Standard 
Model (SM) of particle physics. All observed particles and forces (apart from grav- 
ity) are encompassed by this theory. The SM is a relativistic quantum field theory 
which organizes elementary particles and the forces which act upon them according 
to the principles of gauge symmetry p5fl. Despite its success, though, the Standard 
Model has not furnished a satisfyingly complete explanation of the realm of elemen- 
tary particles. There are many ad hoc parameters, and difficulties arise when the 
theory is extrapolated to energies much beyond a TeV (in the sense explained below) 
without some modification. The Standard Model also affords no link between the 
gauge forces and the force of gravity. 



Supersymmetry (SUSY) |82|, 81, 65 1 is a proposed new symmetry of particle 



interactions which generalizes the symmetries of space-time. The usual space-time 
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symmetries are encoded in a set of commutation relations for the generators of 
translation P m and rotation M mn ^ 

[M mn , Mpq\ — 9npM m q Qmp^nq 9nqM m p + 9mq-^np, 

(1.1) 

[Miwii Pq] — 9nqPm 9mqPni [-fmi Pn] — 0, 

with g mn the metric tensor. These commutation relations constitute a Lie alge- 
bra for the operators, called the Poincare algebra. With reasonable assumptions, 



Coleman and Mandula[45| showed that any Lie group that contains both the Lie 
group underlying the Poincare algebra (the Poincare group) and an internal symme- 
try group must be just a direct product of the Poincare group and that symmetry 
group. There is thus no non-trivial extension of the Poincare group, and one can say 
that the Poincare algebra is the most general Lie algebra of space-time symmetries. 

The notion of a Lie algebra can be extended to that of a graded Lie algebra by 
admitting anticommutators in addition to commutators. The SUSY algebra comes 



from adding so-called "super symmetry" generators Q a to the Poincare algebra (1.1) 
satisfying the (anti)commutation relations 

{Q a , Q b } = i.T) ab Pm, [Q\ Q b ) = o 

[P m , Q a ] = 0, [M mn , Q a ] = \ (a mn Q) a , 
where 7 and a are the usual elements of the Dirac calculus. Haag, Lopuszanski and 
SohniusfMl extended the work of Coleman and Mandula to show that the SUSY 



algebra of equations (1.1) and ([L4 ) constitutes the most general graded Lie algebra 
of space-time symmetries. 

The SUSY generators in ( |1.2| ) are spinorial. Spinors, though non-classical, 
are indispensable in field theories of matter as we know it. They furnish the repre- 
sentation for electrons (and all matter particles), but under the algebra ( |1.1| ) they 
are contingent objects in space-time. SUSY builds spinors into the structure of 



space-time itself and so confers on them an ontological status more befitting their 
fundamental role. Metaphysically, this is tidier. 



In the usual SUSY scheme the transformations of equation |1.2| are global. 
If we promote these to a local symmetry then remarkably we derive the general 
coordinate transformations of general relativity, and find spin-2 gravitons in the 



theory [49] [60]. Local supersymmetry necessarily leads to quantum gravity. No one 
has yet been able to write down a consistent field theory of quantum gravity to 
capitalize on this fact, but it is intriguing and suggests that SUSY may have a 
basic role to play in the unification of the gauge forces with gravity. The glamorous 
"superstring" theory, which is a finite theory of quantum gravity, is founded on 
space-time supersymmetry. 

The notion of unification is on a par with symmetry as a fundamental organi- 
zational principle in physics, as exemplified by Newton's identification of terrestrial 
and planetary gravity, Maxwell's synthesis of electricity and magnetism into a single 
field, or Einstein's equivalence of gravitational and inertial mass. In particle physics 
we can point to the unification^] of the electromagnetic and weak forces in Glashow, 
Salam and Weinberg's electroweak theory, which is a cornerstone of the Standard 
Model. Physicists now anticipate a Grand Unified Theory (GUT) which can unify 
the strong force with the electroweak theory, and a "Theory of Everything" (TOE) 
which will unite these with gravity. 

In the simplest GUT first proposed by Georgi and Glashow (follow- 
ing some work of Pati and Salam[84]), the gauge groups are unified as SU(3)qx 

x Many do not count GSW theory as a proper unification, since one needs as 
many parameters with the theory as without; it is, however, a satisfying conceptual 
unification that puts the electromagnetic and weak forces on a common footing and 
serves as a model for GUT unification. 
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Figure 1.1: Gauge coupling unification in (a) the Standard Model and (b) the MSSM with 



sparticle threshold at 1 TeV. Figure reproduced from 32 



SU(2)l x U(l)y C SU(5) . The running couplings of the three low-energy gauge 
groups are predicted to converge at a high scale where the full SU(5) symmetry is 
restored. Precision measurements of the gauge couplings at Mz by the CERN e + e~ 
collider LEP are consistent with supersymmetric SU(5)|22|, |32|] (because of quantum 
corrections due to SUSY particles lying in mass between the weak and GUT scales), 
but not with minimal nonsupersymmetric grand unification. 

This fact is suggestive. But whether or not SU(5) is the unification group, and 
even if nature "unifies" without grand unification, as in string models, SUSY theories 
generically fix a problem common to nonsupersymmetric theories. This flaw, the 
hierarchy problem |92|], comes from quadratically divergent quantum corrections to 
masses of scalar particles such as the SM Higgs (since there is no symmetry in the 
SM to protect the scalar masses the way chiral and gauge symmetries, respectively, 
protect fermion and vector boson masses). This divergent contribution is due to 
the diagrams shown in Fig. |1.2| . If there were no new physics whatsoever beyond 
the SM, then this divergence would be uncontrolled. So the SM must break down 



at some scale, where new physics intervenes. Very general considerations suggest, 
for instance, that field theories in which gravitation is not incorporated cannot be 
extrapolated beyond the Planck scale (~ 10 19 GeV), where gravitational interactions 
are of the same strength as gauge interactions. New physics must enter here, or 
perhaps at the somewhat lower gauge coupling unification scale ~ 10 16 GeV. 



Perturbative unitarity arguments |50|, 76] require the Higgs mass to be less 
than a few hundred GeV. Keeping the Higgs mass at this scale while integrating 
the quadratic divergence up to, say, A ~ 10 16 GeV would require introducing a 
counterterm into the Lagrangian that must be fine-tuned to an accuracy of (100 
GeV/A) 2 ~ 10 -26 . While not a logical impossibility, this is repugnant to many 
physicists. In order to avoid an unnatural fine-tuning, the scale of new physics can- 
not be much larger than A ~ 1 TeV. If supersymmetric particles have masses at this 
scale, then they can do the job since by the nonrenormalization theorem of Gris- 
aru, Rocek and Siegel]62|], the scalar masses in SUSY theories are not quadratically 
divergent even if SUSY is spontaneously or softly broken^]. The quadratically diver- 
gent contributions from scalars in loops are exactly cancelled by contributions from 



the scalars' superpartners above the SUSY restoration scale, illustrated in Fig. 1.2 . 
Supersymmetry is the only symmetry known which offers generic protection to the 
scalar masses from the quadratic instability to radiative corrections. 

Of course, supersymmetry is not the only extension to the SM which has 



been proposed. Popular alternatives in the literature include technicolor [92] and 

compositeness(85|]. These are attractive theories but they seem to be ruled out in 

2 "Soft" SUSY breaking terms are just those that do not lead to the reappearance 
of quadratic divergences in the scalar masses; these include explicit mass terms 
for scalars (e.g., to lift the quark/squark degeneracy), trilinear scalar interactions 
(important in the stop sector), and others p^j. 
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Figure 1.2: Diagrams relevant to the SUSY nonrenormalization theorem. The scalar loops 
diverge quadratically, as do the fermion loops. Under supersymmetry these contributions 
cancel, and the result diverges only logarithmically. 

their simplest forms. They are strongly coupled theories, like QCD, and are often 
not amenable to calculation^. SUSY on the other hand is a weakly coupled theory 
whose couplings mainly derive from the SM gauge couplings in a canonical way. 
SUSY diagrams are calculable in perturbation theory, and in fact a wide range of 



computer tools exists for this purpose (27|, [79|, 78]. 

That supersymmetry addresses such a wealth of physics issues with such 
economy brings to mind Hertz' comment on Maxwell's equations quoted by Freeman 



Dyson M (from |80|): 



One cannot escape the feeling that these mathematical formulae have 
an independent existence and an intelligence of their own, that they are 
wiser than we are, wiser even than their discoverers, that we get more 
out of them than was originally put into them. 

1.2 The minimal supersymmetric standard model 

An immediate consequence of the spinorial nature of the supersymmetry generator 
is that SUSY transforms bosons and fermions into each other. A supersymmetric 



3 the unreasonable ineffectiveness of mathmatics in the natural sciences, to para- 



phrase Eugene Wigner[97]. 



7 



left-handed 
matter 



right-handed 
matter 



Higgs 
sector 



electroweak 
gauge bosons 



gluon 



spin 



UL 



?L 



UR,d R ,e R 



H° d 



H d J ' ( Hi ) 



spin 1/2 
UL \ ( VeL 



UR,d R ,e R 



H u ,H d ,B, W 



Zi_ (i = 1,2,3,4), 
Wf (j = l,2) 



5 



spin 1 




5 



Table 1.1: MSSM particle content. Only the first generation of matter particles is 
shown; the second and third generations are replicas of this. 

theory must contain for each fermion a boson with identical quantum numbers, and 
vice versa. Spin-| matter fermions get spin-0 partners, and the spin-0 Higgs bosons 
and spin-1 gauge bosons get spin-^ partners. The minimal extension of the SM 



needed to do the job has the particle content listed in Table 1.1. This generic SUSY 
model is called the Minimal Supersymmetric Standard Model (MSSM). Supersym- 
metric particles are marked with tildes. By convention, fermionic superpartners of 
SM bosons take the suffix "ino," and scalar counterparts of SM fermions have the 
prefix "s" (short for "SUSY"). The term "sparticle" refers to any of these super- 
symmetric particles. 

Matter fermions are represented by Dirac spinors and have four degrees of 
freedom, while (complex) scalars have only two, so each fermion / is associated with 
a pair of scalars. Conventionally these are chosen to be the sparticles corresponding 
to the fermion's two helicity states, and are written as /l and f R . (This notation 
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is perhaps confusing since scalars of course have no chirality; it is still true, though, 
that only fi couples to the W (as also only ft)-) The particle studied in this thesis 
is the lighter supersymmetric partner of the top quark (a linear combination of ti, 
and tfi), variously called a top squark, scalar top, or stop. For brevity, we use "stop" 
in the present work. 



The only new non-super symmetric content in Table 1.1 is an extra Higgs 
doublet, placing the MSSM in the class of Two Higgs Doublet Models (2HDMs). 
In the Standard Model, three of the Higgs doublet's (real) degrees of freedom go to 
supply longitudinal polarization for the electroweak gauge bosons, leaving a single 
neutral scalar Higgs. In 2HDMs similarly 8 — 3 = 5 states remain for physical 
particles: two CP even Higgs bosons, a CP odd Higgs boson, and a pair of charged 
Higgs. 

This extra Higgs doublet is a consequence of consistently applying supersym- 
metric field theory. The technique involved — the so-called superfield formalism — is 



too lengthy to go into here, but many standard treatments exist (see [70j and ref- 
erences therein). One writes down the most general supersymmetric Lagrangian 
including soft SUSY breaking terms (which lift the mass degeneracy of superpart- 
ners) consistent with Poincare and gauge invariance. With this treatment it is easy 
to see that SUSY requires separate Higgs fields to couple to up-type and down- 
type quarks. An immediate consequence of this fact is the introduction of two new 
fundamental parameters. One is the ratio of the VEVs for the two Higgs scalars, 
tan/3 = Vw/v^. The other is a higgsino mixing parameter [i which parameterizes 
the coupling ^lH u H,i; [x is required to be non-zero to prevent problems due to a light 
higgsino. 
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The fermionic partners of the SM photon (photino), Z boson (zino) and two 
neutral Higgses (higgsinos) all share the same SU(3)xU(l) e m quantum numbers 
and will, therefore, mix if electroweak symmetry is broken. The resulting states, 
called neutralinos, are denoted Z\, Z2, Z3 and Z4 in order of increasing massQ. For 
many phenomenological considerations, it is significant whether the Z\ has a pre- 
dominantly gaugino or higgsino character. Since we typically set supersymmetric 
branching fractions to 100% for the decays we study here, this distinction is largely 
unimportant for our work. As with the neutralinos, the superpartners of the W and 
the charged Higgs also mix to form two chargino states, the lighter W\ and the heav- 
ier W2- The fermionic content that the MSSM adds to the SM (gauginos+higgsinos) 
is free of chiral anomalies. 

General techniques for the construction of supersymmetric particle Lagrangians 



are discussed in [93]. In their most general form, these SUSY Lagrangians admit 
gauge-invariant renormalizable lepton- and baryon-number violating interactions, in 
contrast to the Standard Model. The completely unconstrained MSSM scheme is 
minimal with respect to particle content, but maximal in terms of parameters — 106 



new parameters are required in addition to the SM's 18 [52]. These new parameters 
are strongly constrained by current data on CP violation, flavor-changing neutral 
currents, and lepton flavor non-conservation. Setting these possible phases are all 
to zero, the resulting theory has 30 new parameters. The MSSM referred to in this 

thesis is this 30+18 parameter model. 

4 There are two systems of notation for particles in the electroweak gaug- 
ino/higgsino sector. In this thesis we call the neutral members of this group Z{ 
(i = 1,2,3,4) and the charged members (j = 1,2) (although we usually sup- 
press the charge sign and just write Wj) out of prejudice for the gaugino character 
of Zi t 2 and W\ expected in the mSUGRA models introduced below. In other works 
these neutral and charged "inos" are written as Xi an d xt > respectively. 
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B and L conservation leads to a further discrete symmetry of the renor- 
malizable terms in the Lagrangian, imparity invariance, with quantum number 
R = (—i~) 3B+L+2S ( w here S = spin). All SM particles (plus the second Higgs 
doublet) are i?-even, and all superpartners are -R-odd. If i?-parity is a good sym- 
metry then sparticles must be pair-produced by collision of SM particles, and a 
decaying sparticle must have an odd number of sparticles among its decay products. 
Therefore there must be a lightest supersymmetric particle (LSP) which is stable. 
Limits on long-lived colored or charged particles require the LSP to be weakly in- 



teracting and electrically neutral] 54]. The Z\ is generally taken to be the LSP. (The 
only other LSP candidate is the scalar neutrino u, but the sneutrino is disfavored 
if the LSP is also to constitute the cold dark matter presumed to lie in the galactic 
halo|72|j.) In SUSY events at particle colliders, this LSP will escape the detector 
without interacting, leading to the classic SUSY signature of missing transverse 
energy # T . 

If SUSY is to solve the hierarchy problem, then the masses of the super- 
partners of the MSSM cannot be much greater than 1 TeV, and so particle colliders 
currently in operation (especially the Fermilab Tevatron) or under construction (the 
CERN LHC) have a good chance of finding some evidence for their existence. This 
accessible scale of new particle states is small compared to the scale imagined for the 
mechanisms underlying supersymmetry (from 10 10 to 10 19 GeV), and these theories 
are generically called "low-energy supersymmetry." 

1.3 Constrained SUSY models 

The MSSM is completely agnostic about physics above the TeV scale of new particle 
states. The many parameters mentioned in the previous section are so-called "soft" 
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SUSY breaking parameters, and are written down from very general considerations 
of the supersymmetric Lagrangian. The parameter space of the general MSSM is too 
large for viable phenomenological studies, and one therefore turns to more restricted 
models. Further constraints among these parameters require a theory of physics at 
energies beyond 1 TeV. There being a paucity of data on physics at such scales, no 
complete and compelling supersymmetric theory has emerged as yet. 

Constraints can be found by considering the nature of supersymmetry break- 
ing. SUSY must of course be broken since no sparticles degenerate in mass with 
their ordinary partners (e.g., scalar electrons) have been observed. SUSY breaking 
relying only on the MSSM content of the theory (or even on just the TeV scale 
with more particles) seems to be infeasible, so models invoke a "hidden" sector of 
particles which have no SM gauge interactions. In supergravity models, SUSY is 
broken in this hidden sector, which communicates with the visible sector through 
gravitational interactions alone. One obtains a simple structure which one assumes 
is valid at a high energy scale (Mqut or -^planck) an< ^ uses ^normalization group 
equations (RGE) to calculate low energy MSSM parameters. The high energy pa- 
rameters are usually universal (for instance with gravity-mediated SUSY breaking), 
and a GUT ansatz affords further simplification by relating gaugino masses. 

The minimal supergravity (mSUGRA) models widely used in phenomenolog- 



ical studies are prepared according to this prescription [56|. mSUGRA depends on 
only four continuous parameters and a sign in addition to the 18 parameters of the 
SM. These are generally taken to be a universal scalar mass mo, a universal gaugino 
mass rrii/2, a universal trilinear interaction Aq (sfermion-Higgs-sfermion) , the Higgs 
VEV ratio tan /3, and the sign of the Higgsino mixing parameter fi: 

mo, m 1/2 , A , tan/?, sign Qu) (1.3) 
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all input at the GUT scale except for tan/3, which is a weak-scale quantity. The 
Higgs bosons, which are part of chiral supermultiplets, share the common scalar 
mass mo at the high scale, and are renormalized due to gauge interactions like the 
doublet sleptons. The Higgs doublet that couples to the top family gets very large 
negative contributions to its squared mass rn 2 H due to the large top Yukawa coupling 
(discussed below) which can easily be driven to negative values leading automatically 
and naturally to radiative electroweak symmetry breaking (EWSB).0 In mSUGRA 
the magnitude |/x| is fixed by requiring EWSB to reproduce the measured Z boson 
mass. 

mSUGRA has many distinctive features, and its phenomenology has been 
studied extensively. We collect a few facts here which will be useful in our work. 
Because the high-scale mass parameters are universal and the RGE depend mainly 
on the known gauge couplings, one gets generic mass relationships such as 

m~ = Trig + m 2 q + (5 — $)m\i 2 + -D-terms (1.4) 

where the D-term scale is M|/2. The squark, slepton and gluino masses obey 

m | = mj + (0.7 - 0.8) mf (1.5) 

(ignoring quark and lepton masses). The lighter neutralinos and lighter chargino 
are mostly gaugino (rather than higgsino) and one generally has 

2m~ « « m-g/3. (1.6) 

The relations (|0]) hold because gaugino masses are assumed to unify at the GUT 
scale. 

5 The electroweak symmetry breaking is automatic and natural, but not miracu- 
lous; one still has to input the weak scale by hand into the high energy parameters 
of the theory. 
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mSUGRA is just one of a large and growing set of models intended to capture 
physics at the high scale. There are theories that are more constrained, such as cer- 
tain superstring-inspired models, in which only one or two free parameters appear. 
There has also been recent interest in models which vary the different mSUGRA 
assumptions on universality at the GUT scale. For instance, this may happen if 
masses unify at M~p\anck anc ^ become non-universal as they run down to ^QUT> 
or if the unifying group is larger than SU(5) and admits additional mass contri- 
butions. String models, hypercolor models, and models with novel SUSY breaking 



mechanisms can also exhibit non-universality. Reference [23] surveys the gross phe- 
nomenology of these models. Because the light stop sector depends on relatively few 
MSSM parameters, our results are not particularly sensitive to the actual physics 
at the high scale. Our results, as we will see, will broadly depend only on how the 
stop is assumed to decay. 
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Chapter 2 



The Top Squark 



2.1 The light stop 

The top quark has recently been detected and its mass measured in the lepton 
+ jets||, i 1, dilepton|5|, g] and all-jetsQ channels. A concise summary of 
these results is presented in [|7l| which derives a combined value of rat = 174.3 ± 

3.2 (statistical) ±4.0 (systematic) GeV. For the purpose of our calculations, we adopt 
m t = 175 GeV. 

With such an extraordinary mass, the top family Yukawa interactions are 
large and comparable to the electroweak gauge interactions. In renormalizing squark 
masses from the unification scale down to the physical mass scale, this top Yukawa 
coupling depresses the diagonal masses of the left and right stop states ii and Ir 
relative to the masses of the generation 1 and 2 squarks and also allow substan- 



tial mixing between the two states [58, 69, 68]. The stop masses come from Higgs 
superpotential terms 
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(h^ and h° d are the up- and down-type Higgs superfields, t and T c are left and right 
top superfields) and soft SUSY breaking terms 

£ soft SUSY breaking 3 A tfthH^ R - ml$ L i L - m 2 lR $ R t R . (2.2) 
The stop masses are (assuming universality as in mSUGRA) parameterized 



as 



GUT-scale mass gauge RGE Yukawa RGE EWSB mass 

m\ L ml + Fm\ /2 - 2 (/ 2 A t + / 2 A 6 ) m 2 + m 2 t 

D term 



-(-cos2/3)Q-^sin 2 ^) Mf (2-3) 
= ml L +m 2 -2(f 2 A t + f 2 A b )ml 

m \ R = m L + m t - 4 a 2a *^o- 

where the I s * and 2 nc ^ generation left- and right-handed squark masses at the weak 
scale are 

m l L = m l + Fm i/2 ~ (~ cos sin2 °w) M|, 



m\ R = ml + Fm 2 /2 - (- cos 2/3) (§ sin 2 6» w ) M 2 



(2.4) 

z- 



At the GUT scale the stops get contributions from the universal scalar mass tuq. 
Gauge interactions in the RGE (predominantly QCD) add the term proportional to 
the gaugino mass, where the coefficient F has a typical value of 5 — 6. The top and 
bottom Yukawa couplings are 

ft = 9 2 m 2 t/ (2M 2 , sin 2 /?) « 1.01/ sin 2 /?, 

(2.5) 

fl = g 2 m 2 b / (2M^ cos 2 0) 0.00082/ cos 2 

with g the weak coupling. The quantities A t f, are calculated from the RGE and 
take values on the order of 0.1. For small tan/3 the top Yukawa ft dominates and 
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m| < mi h since Ir gets twice as much depression as Il\ these two masses approach 
each other as tan/3 — > ratjra^. The usual EWSB mechanism adds the top- flavor 
mass term mf. The factor cos 2/3 = (l — tan 2 /3) /(l + tan 2 /?) in the L>-term is 
negative for tan/3 > 1, so the Z)-terms depress the diagonal masses. 

The Yukawa interaction also mixes tj_, and tR and the stop mass matrix from 



interactions (£DJ) and ( gg ) isg§ 



T7T-7 



(2.6) 



with at = — /z cot /3. Diagonalizing this gives mass eigenvalues 

n 2 



11,2 



m~ + m 



'2 2 \ 2 

m~ — m| \ 



\ 



3//A t + /^A 6 



(2.7) 



i?A t m 2 



1 - - sin 2 9 W 



cos 2/3 



+ (A t — n cot /3) 2 m 2 



and the physically propagating eigenstates 



cos at — sin ft 
sin t cos 4 



t L 

tR 



rotated by the stop mixing angle 

tan6*t = 



o 

ra- 
th 



2 

ti 



2.J 



mt (-4* - M cot /3) 

In the absence of mixing, Ot — > vr/2 and ti = tR in models with a universal mass at 
the high scale. 

For moderate tan/3 the stop mass splitting is dominated by the off-diagonal 
terms in the mass matrix. The unification scale trilinear coupling parameter Aq 
can be freely adjusted, and light stops can usually be arranged irrespective of the 
value of the other universal mass parameters tuq and raj^; such parameter sets can 
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generally be made consistent with experimental constraints (particularly those on 
Mz and m#)[|5l|]. In fact, t\ masses can be driven low enough that model builders 
have to be careful to insure that the stop does not become the LSP or, worse, 
tachyonic leading to dangerous color-breaking minima. 

2.1.1 Theoretical prejudices 

We have seen that it is possible for the top squark to be light. There are also 
arguments that directly favor a light t\. Chief among these is electroweak baryoge- 
nesis, a mechanism for generating the observed baryon number of the universe at 



the electroweak phase transition JTJ]. The original baryogenesis idea of Sakharov[88] 
involving baryon number violation, C and CP violation and thermal nonequilib- 
rium seems not to work in the Standard Model. The CP-violating phase in the 
Kobayashi-Maskawa matrix is too small to allow enough baryon number production 
in the symmetric phase, and experiment rules out a Higgs scalar light enough to 
generate a strongly first order phase transition which would lock in the generated 
baryons. SUSY, though, has many extra CP-violating phases. And a light (scalar) 
stop strongly coupled to the Higgs through its large Yukawa coupling can provide 
for the first order phase transition while allowing the Higgs to be heavier. De- 
tailed calculations [43, 48] indicate that the stop mass must lie in the range 100 GeV 



< 



~ m t 

find that this range is quite accessible to the Fermilab Tevatron using the searches 
outlined in this thesis. 

Let us collect here a few other theoretical considerations related to the stop. 
First, as mentioned above, the LSP Z\ is a favored candidate for the dark matter 
determined to be present in the halo of our galaxy. In order for these relic neutralinos 



l x < 160 GeV for this mode of electroweak baryogenesis to succeed. We will 
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Figure 2.1: contributions to the decay b — > from the Standard Model (SM), two Higgs 
double models (2HDM), and supersymmetric models (SUSY). 

to have an appropriate density, it is necessary that they cannot have annihilated via 
Z\ + Z\ — > SM particles at too great a rate. If t\ is not much heavier than Z±, then 



it is also possible to spoil the relic density through co-annihilation 1 38]. The authors 



of [38] have studied this process, and concluded that maintaining the favored relic 
density 0.1 < £1^ < 0.2 requires a mass gap mr — m~ > 11 to 33 GeV. Excluding 
the small mass gap is good news for collider experiments, where small — m~ 
means poor detection efficiency. 

There are also constraints on the stop sector coming from precision elec- 



troweak measurements. One of these is the p parameter [55], which is a measure 
of the ratio of charged current to neutral current interaction strength. The close 
agreement between LEP II experimental determinations and the Standard Model 
Ap = means that ti, which couples to the W, must be quite heavy. A recent 



study H finds that in mSUGRA one needs m~ th > 275 (310) GeV for p > (/x < 0). 
p parameter limits on the W\ and Z\ masses are weaker than the current direct 
experimental limits. 

In 2HDM models with a light charged Higgs the rare decay b — > s~f may 
require a light stop as well. This decay has been measured[21] to have a rate close 



to its SM predicted value. A light H enhances the decay by allowing a charged 
Higgs diagram analogous to the SM W~ diagram, as in Fig. ^O], while a SUSY 
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diagram with a chargino and light stop in the loop decreases the rate and can bring 



it back into line with the SM rate[31, 83]. This is generally not a problem in models, 
as the charged Higgs is typically heavy. However, the p parameter data allow an 
mSUGRA H~ as light as 140 GeV (fi < 0) g|. 

2.2 Stop production and decay 

Stops are color triplets, as are top quarks, and at a hadron collider such as the 
Fermilab Tevatron will be strongly produced in t*t% pairs (to conserve R-parity) 
through gluon fusion and qq annihilation. The tree level production cross-section 



(Tj,^ depends only on the mass m^. Figure 2.2 shows this tree level o^-*^ as a 
function of for the Fermilab Tevatron pp collider at a center of mass energy 
y/s = 2.0 TeV. These cross-sections were produced by the program ISAJET (see 
below) using the CTEQ2L[] parton distribution functions |75[| . For comparison, the top 
quark pair production cross section for 2.0 TeV pp — > ti is a t i = 6100 fb when 
calculated by the same means. 

The Next to Leading Order (NLO) corrections to the tree-level o^*^ have 
been calculated]^]. The authors of [|36| used mSUGRA-inspired parameters to 
calculate higher order corrections for the Tevatron pp collider at ^/s = 1.8 TeV. The 



1 We have checked that our results are not affected by using more modern parton 
distribution functions than the CTEQ2L of our original analysis; this is as expected 
since the relatively heavy stop pairs are produced at high x. 



m h (GeV) 


K 


70 


1.41 


110 


1.30 


150 


1.19 


190 


1.11 



Table 2.1: K factors for Tevatron stop pair production. 
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Figure 2.2: Stop pair production cross-section at the Tevatron. 



corrections are positive and substantial. Table 2.1 lists the K factors they obtained. 
We use the lowest order cross-sections for both signal and background in the present 
work. 

Current experimental mass limits for SUSY particles of interest here are listed 



in Table 2.2. These are generalized lower bounds which we will use for reference in 



this thesis. Many of these limits are correlated, and in some cases higher masses 



particle mass (GeV) 



notes 



MM 
1,0 

large mo,® © M\ 
large mo,(l£, 11] 

in 

E3 



a 
h 

Wi 
Zi 
e 
H 



220 

173 

90 

93 

32 

92 

108 



Table 2.2: Current mass bounds for particles relevant to this study. 
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Figure 2.3: Decays of the light i\ studied in this thesis. 

have been excluded for particular combinations of other SUSY parameters. It is 
also possible to find corners of parameter space where these bounds may be evaded. 
See the indicated references for a fuller discussion. The t\, W\ and Z\ experimental 
determinations are discussed briefly in Section 1.2.3. 

Given these masses, the possible light t\ decays in SUSY models without 



exotic particle content are [24] the 2-body tree-level modes t\ — ► tZ\ and t\ — > bW%, 



the 2-body loop mode t\ — ► cZ\, the 3-body tree- level modes t\ — > bWZ\, t\ 



and t\ — > biP, and the four-body decay t\ — > bf f Z\. Figure 2.S diagrams the decays 



analyzed at length in this thesis, and Fig. displays the others. Table O presents 
a summary of these modes. 

We could write decays analogous to those of the preceding paragraph with 
a gluino g replacing the LSP Z\, such as t\ — > eg. The current gluino mass 
bound (see Table |2.2|) would require a rather heavy stop to open this mode and 
overcome phase-space suppression enough to compete with the (loop-suppressed) 
t\ — > cZi, especially under the gaugino unification relationship rrig ~ 6m~ . Even 
if the reaction did have an appreciable branching fraction, it would be attended 
by the cascade g — > qqZ\ and so would be similar the decay t% — ► cZ\ (which we 
study at length) but with more jets and a softer LSP. The other possible gluino 
modes are even more suppressed — t\ — * bWg has too heavy a final state and 
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Figure 2.4: Decays of the t\ not studied in this thesis. 

t\ — > bff'g has too little phase space. We do not consider gluino decays fur- 
ther^ 

When > mt + the direct 2-body decay ti — ► tZ\ is open, but it will 
be strongly suppressed by phase space unless the mass gap is large. This requires 
mi 3> 210 GeV. From Fig. |2.2| the production cross-section for such stops is less 
than 150 fb, and at l/40th the ti production cross-section the likelihood of observing 
this mode is remote. We do not consider it further. We also mention in passing that 



the supersymmetric top decay t—*t\Z\ can open up [89] if the stop and neutralino 
are light enough. See section 1.2.3 for this. 

If > mft + , the 2-body tree-level decay t\ — ► bW\ will dominate. 
The chargino from t\ — > bW\ rapidly decays to a neutralino and fermion pair, 
W\ — > f'fZi, so that t\ pair production results in t\t\ — > bb + 4 fermions + Z\Z\. 
(We do not consider models with a slepton light enough for W\ — > £u, as explained 

below.) This top squark event topology is similar to the Standard Model ti pattern 

2 Models with light gluinosj25|, Q which could have evaded detection have been 
proposed in the literature. We do not consider such models here. 
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ti — > 66 + 4 fermions, the main difference being the presence of the Z\Z\ pair. We 
will find that for the range of t\ masses accessible with the luminosity upgrades at 
the Tevatron a real chargino cannot decay into a real W boson, which provides for 
another difference from ti events. (The chargino will decay to a real W only when 
niyy — m~ > mw- Given the measured W mass and the experimental limit on 
the neutralino mass, this requires > 115 GeV, rising to > 160 GeV as 

— ► 2m^ (see eqn. |Q| ) and we use the limit. As we will see in Chapter 
H such heavy charginos are not accessible at the Tevatron in i\ — ► bW\ decays.) 

Phenomenologically, t\t\ events for the t\ — > bW\ mode are categorized 
according to the leptonic or hadronic nature of the chargino decays just as ti events 
are categorized according to the nature of their W decays. There are (1) dilepton 
events t-fti — ► bb£ + £'~ vvZ\Z\ whose signature is (6)-jets + dilepton + |? T , (2) 
1-lepton events t\i\ — » bbq' 'q£vZ\Z\ with signature (6)-jets + lepton + ^ T , and 
(3) all-jet events t\t\ — * bbq' qq" q'" Z\Z\ with (6)-jets + The all-jets channel 
suffers from large QCD backgrounds; the search is difficult in this channel and 
we do not consider it further in this study. The leptonic channels show the same 
topologies as the corresponding tt channels, and thus the important backgrounds for 
these channels are those usually identified in ti studies, along with ti itself. These 
channels are considered in detail in Chapter || 

When stop decays via the chargino t\ — * bW\ then one also has to consider 
the cascade decay W\ —*ff'Z\. We are particularly interested in the branching 
fraction to leptons, B (W\ — > leptons + Zij , since our detection schemes look for 
these leptons. When the decay proceeds predominantly via a virtual W, as in 
Fig. |2.5| a, then the chargino branching fractions will be those of the mediating W, 
B(Wi -► e) = B(Wi ->(!,) = B(Wi -► t) = B(W -»■ e) » 1/9. But if some sfermion 
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f 

(a) (b) 



Figure 2.5: Tree-level decays of the chargino. The diagram at left is the W^-mediated 
decay, which will lead to VK-like fermionic branching fractions if it dominates. The sfermion 
diagram on the right will enhance the branching fraction to / / when / is light. 



/ is light enough for the diagram in Fig. |2.5| b to compete, then the branching 
fraction to the corresponding fermion and its partner W\ —tff'Zi can be enhanced. 
If sleptons are lighter than squarks (but still heavier than W\) then we can get the 
branching fraction B (Wi leptons + ZA » 1/9. In mSUGRA (and in many other 
models) we can easily have iriq 3> from eqn. ( |1.5[ ), and one typically finds an 
enhancement. This almost always works in our favor, as discussed below in Chapter 
& The /'"-mediated decay can also be important when the WW\Z\ coupling is 
dynamically suppressed. 

In case < m& + , the stop must go via the 2-body loop decay or the 
3- or 4-body decays. In mSUGRA > and experimentally > 80 GeV 



so we will not consider t\ — > btv or t\ — * biv here; they are treated in [47, 51]. 

When mi x < nit, + mw + m z l > ^ ne remaining 3-body decay is closed, and t\ 
must decay as t\ — > cZ\ or t\ — > bf f Z\. The rate for the loop decay t\ — > cZ\ was 



estimated in pq], and earlier workpl, 28| had focused on the this decay. Recent 



calculations [39] indicate a substantial branching fraction for the 4-body decay for 



realistic regions of the MSSM parameter space. The event topologies for t\ 
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cZ\ 

U V V 1 


discussed fully in Chapter ||| 

H 1 CQf^H Ti ill "\r m 1 , M PI T^1~ pr 1^1 
<_LlDi^ U-OOCVJ. 1 Lll-Lj HI V^lldJJ IC1 lUI 


bWZ 1 
b£u, btv 
bff% 


competes with ii — > cZi when 3> 165 GeV; see 86, jS7f 
requires light slepton; not studied here; see [^7], [5l|| 
may compete with t\ — ► cZ\\ see [|39|] 


eg 
tz x 


needs mr S> 180 GeV and loop-suppressed 
needs m~ ti > 210 GeV 



Table 2.3: Light stop decay modes. 
bff'Zi, but not kinematics, would be similar to those for the t% — > 6Wi decay. It is 



interesting to note that Ref. 3E] finds a gradual transition of dominance between 
the two decay modes, so that there are broad regions where both decays occur 
appreciably. This of course would greatly complicate the search strategy and make 
detection more difficult. In this work we only consider the loop decay titi — ► 
ccZ±Zi, whose topology is c-jets + $ T . The classic $ T + jets signature suffers from 
big backgrounds from W and Z events with jets from QCD radiation and^ T supplied 
by neutrinos. We will see below that charm tagging is necessary to develop good 
detection schemes. Chapter || is devoted to searches in this channel. 

In case the t\-Z\ mass difference is great enough to open up the 3-body 
decay, then t\ — > cZ\ and t\ — ► bW Z\ will compete — the 2-body decay is loop 
suppressed and the 3-body decay is suppressed by phase space |8^, p7| . Note that 
for the common ansatz = 2m. ~ discussed earlier, the 3-body decay only opens 
up for mi x > 165 GeV (with > 160 GeV), and even then is strongly suppressed 
by phase space so this mode can only be important for 3> 165 GeV. This 
region is difficult for either the t\ — > cZ\ or t\ — > bW\ search strategies we develop 
below for 2 fb _1 integrated luminosity at the Tevatron, and a direct search for 
t\ — > bW Z\ would be even more problematical due to the presence of on-shell Ws 
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(whose rejection is key to suppressing SM backgrounds). Our preliminary analysis 
has indicated that this decay mode is inaccessible to the Tevatron experiments, and 
other investigators have since found the same result. We do not pursue the 
3-body decay further. 

Figure [2~^ shows the regions where these various decays are expected to occur. 
The neutralino decay mode t\ — ► cZ\ and the chargino mode t\ — ► bW\ are overlaid 
in this figure. The two axes will coincide if = 2m^. The upper left region 
marked < is ruled out because Z\ is the LSP by assumption. 

We restrict our attention here to physics at the Fermilab Tevatron pp collider 
(described in Section 1.3.1). There also exists a rich stop phenomenology at e + e~ 
colliders (see |35|] and references therein), and proposed fi + colliders pi]. We do 
not examine .R-parity violating processes for the Tevatron, which are discussed in 



IH and PI 



A note on the top sector 

Since we are interested in assessing the reach for the SUSY signal over Standard 
Model backgrounds, we assume that the top always decays via t — > bW and neglect 
SUSY decays of tops such as t —* t\Z\ f67|, [T?], or t — > t\g\9C}. If these decays 



are present to any extent, they could potentially add to the signal as the stops 
cascade via the decays discussed above. That is, a tt event where one of the tops 
decays supersymmetrically will be less likely to fail the cuts designed to reject the tt 
background, and so will improve the yield. Our assumption leads to a conservative 
background estimate, and if anything will underestimate the discovery potential for 
stops. 
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Figure 2.6: Map of decay modes and experimental results for the light stop. The vertical 
axis on the left is neutralino mass, and that on the right is chargino mass. In case rn^ — 
2m~ , these axes coincide as illustrated. The heavy dotted line is the CDF 95% exclusion 

contour for the mode t\ — > cZ\. The dashed line shows the LEP II 95% confidence limit for 
neutralino mass and stop mass given the decay i\ — > cZ\. The dot-dashed line is the LEP 
II chargino mass limit. 
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2.3 Current experimental results 

There are good experimental limits on the stop mass from both the CERN LEP II 



e + e~ collider and the Fermilab Tevatron pp collider. LEP II results) 1C, 17, 29] at a 
collision energy ^/s = 189 GeV set a 95% confidence limit exclusion on stops lighter 
than 90.3 GeV for the decay mode t\ — > cZ\, assuming the stop is at least 10 GeV 
heavier than the neutralino. This limit weakens by a few GeV when 9t = 0.98 and 
t\ decouples from the Z. The LEP experiments do not report a stop mass limit for 
the decay t\ — ► bW\, but their chargino mass limit (see below) implies > 98 
GeV if stops decay via this mode. 

At the Tevatron, the DO collaboration H investigated the t\ — * cZ\ mode 



in the $ T + jets channel (see Fig. 2.6 and Section 2.2) and excluded the region 
bounded by m ix < 100 GeV, m~ < 40 GeV and m~ ti - to- > 40 GeV. The CDF 
collaboration |l^] got a bigger excluded region by analyzing 88 pb _1 of Run I data 
(see Chapter ||) in the-^ T + jet + charm jet channel for this mode and set a 95% CL 
exclusion for a region roughly bounded by < 115 GeV and m~ < 50 GeV. The 
Tevatron collaborations also searched for the ii — * bW\ mode. CDF looked for a^ T 
+ lepton + jets (from b-quarks)[pC|1, and DO looked in the dielectron channel||. Nei- 
ther group was able to improve the existing limits (of > m^. + > 95 GeV) . 

LEP II also provides the best limits on the W\ and Z\ masses. Chargino limits 
are generally function of mo; small mo allows W\ and Z\ to be nearly degenerate 
in mass which seriously compromises the detection efficiency and hence the mass 



reach. The latest published results (for a collision energy yjs < 189 GeV)(Tl|, IS, 3C] 



find > 93 GeV and > 32 GeV for 3> m^r (weakening to m^. > 72 
GeV, > 31 GeV in the worst case). 
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Experiments at LEP II are ongoing, and limits from there are evolving even 
as this thesis is being written. The experiments now have accumulated data at 



\fs = 200 GeV which will published shortly. For instance, ALEPH[53| has just 
released preliminary bounds of rn^ > 100 GeV, > 35 GeV, mg > 92 GeV, 
ms > 85 GeV, m H > 92 GeV, and tan/3 > 1.9. 
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Chapter 3 

Tevatron Simulation 



3.1 Fermilab Tevatron 

The Fermilab Tevatron is currently the world's premier (indeed, only) pp collider. 



It has two major general purpose collider detectors, CDF[|12|, 41] and D0[| 
For reference, we group Tevatron operations into three phases: Run I, Run II and 
Run 11+ . Run I refers to the data already collected during the machine's previous 
operating phase in 1994-1995. Run I obtained an integrated luminosity of almost 
100 pb _1 in each detector at a center-of-mass energy of y/s =1.8 TeV. These data 
have already been analyzed for signatures of the light stop, and the results of these 



analyses are presented below. Table 3.1 lists the parameters for these runs [p4]| . 

The Tevatron is currently undergoing major machine and detector upgrades 
and is scheduled to resume operation in 2001. Chief among the upgrades is the 





(TeV) 


C (cm 2 sec x ) 


fCdt (fb^ 1 ) 


dates 


Run I 


1.8 


< 10 31 


~ 0.1 


1994-1995 


Run II 


2.0 


2 x 10 32 


~ 2 


2001-2003 


Run 11+ 


2.0 


2 - 5 x 10 32 


> 15 


to 2007 



Table 3.1: Tevatron run parameters. 
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new "Main Injector" which will store and pre- accelerate protons and antiprotons for 
injection into the Tevatron ring itself. The Main Injector will also recycle antiprotons 
from the Tevatron resulting in a tenfold increase in luminosity. Together with the 
new Recycler Ring, the peak luminosity is expected to rise by a factor 20 from 
1 x 10 31 cm~ 2 sec to 2 x 10 32 cm _2 sec _1 . Additionally, the beam energy will be 
raised to yfs = 2.0 TeV (which will increase the t*t and ti production cross-sections 
by 40% over those at y/s = 1.8 TeV). The primary data collection period following 
the Main Injector upgrade is referred to as Run II. During Run II we expect that 2 
fb _1 of data will be accumulated at 2.0 TeV after about two years of operation. 

The Tevatron is expected to continue operations for several years beyond the 
nominal 2 fb _1 Run II phase. There are also proposals on the table for further 
major luminosity upgrades. In anticipation of a substantially larger data sample, we 
also make projections for an integrated luminosity of 25 fb _1 at 2.0 TeV, which we 
call Run 11+ . 

3.2 Event Simulation 



Collisions at pp machines are messy affairs [33]. At 2 TeV the proton is a complicated 



assemblage of quark and gluon "partons" interacting via the strong force. In a hard- 
scattering event almost always just one parton from the proton and one parton from 
the antiproton interact. Each scattering parton has a certain probability to be of 
a given type (say gluon, strange quark, or whatever) and to carry a given fraction 
of the hadron's momentum. The set of functions which encode these probabilities 
for given momentum transfer in the hard-scattering sub-process are called parton 
distribution functions. They are fundamental inputs in any model that describes 
large momentum transfer hadron scattering processes.. 
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The partons can be treated as free on their way to the hard-scattering, as are 
their colored reaction products. While inbound they radiate gluons freely, some of 
which are hard enough to be resolved as initial-state radiation (ISR). The partons 
being free at the collision point, the hard scattering reaction itself can be calculated 
in perturbation theory. After the collision the reaction products recoil from each 
other and the ISR. Color confinement causes the colored reaction products, any 
final-state radiation, the ISR, and the remnants of the initial p and p to hadronize 
as they flee the collision point. The usual picture is of color flux lines which stretch 
and break in a process called fragmentation. Each escaping colored particle will 
evolve into a multitude of color singlet hadrons, all moving in rough collinearity to 
the original particle. This concentrated shower of particles is called a "jet." (The 
remnants of the original p and p move down the beam axis to form "beam jets.") 

The hard scattering processes are described by the underlying high energy 
theory. In our case , these are 

99 «*ti 

and are described by the super symmetric model discussed in Chapter 1. The decays 
of the stops into SM quarks are also described by this theory. The processes by 
which the partons fragment and ultimately produce hadronic jets, leptons, photons 
and so on that are detected in the experimental apparatus is independent of the 
high energy theory, and phenomenological models that describe these are encoded 
in event generators that we outline below. 

The (stable) hadrons, charged leptons, or photons produced in the collision 
move away from the collision point and record their passage in the surrounding de- 
tector. The detector is a complicated apparatus designed to characterize the visible 
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reaction products. It has components for tracking charged particles, calorimetry for 
absorbing and measuring the energy carried away from the collision by the hadrons 
as well as electromagnetic particles (mainly electrons and photons), and on the 
outside a system for measuring muon momenta. Modern detectors also include a 
sophisticated system that serves to resolve displaced secondary vertices from the 
decay of heavy flavors. 

The fraction of the hadrons' longitudinal momentum carried by the initial 
hard scattering partons is unknowable, so there is an irreducible uncertainty in 
the longitudinal boost of the center-of-momentum frame for the colliding system. 
However, by forming the vector sum of all the energy deposited in the transverse 
direction, we get an important quantity called missing transverse energy, ^? T . A 
certain amount of^ T is due to jet and lepton mismeasurement from imperfect en- 
ergy resolution, particles going into cracks in the detector, and other "non-physics" 
causes. A large however, generally indicates the production of one or more 
high-energy long-lived weakly-interacting particles that escape the experimental ap- 
paratus without depositing energy. In the Standard Model these would be neutrinos. 
In SUSY searches, large ^ T is the signature of escaping LSPs. Indeed, since a pair 
of LSPs are always produced in a SUSY reaction (for i?-parity a good symmetry), 
large $ T is the hallmark of a supersymmetric reaction. 

3.2.1 Collision simulation 



We use the program IS A JET v7.29[^ with the IS ASUS Yg] extension to simulate 
Tevatron events for both t-fti and background processes. ISAJET does Monte Carlo 
event generation for pp (and other) collisions. It supports a wide variety of parton 
distribution function sets (We use the CTEQ2L set.). It incorporates lowest order 
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perturbative QCD cross-sections, initial- and final-state QCD radiation, indepen- 
dent fragmentation and final-state hadronization, and underlying event effects. It 
includes a phenomenological model tuned to minimum bias and hard scattering 
data for the beam jets. ISASUSY provides supersymmetric particle production and 
branching ratio calculations, and keeps track of all SUSY decay chains. A SUGRA 
package accepts GUT-scale input (Eq. Q1.3D ) and applies the RGE to produce low- 
energy MSSM parameters relevant for phenomenological calculations. 

Because we control the relevant SUSY parameter space directly, we use the 
ISASUSY MSSM technique to set SUSY masses , and by hand. Events 
are generated with the reaction products' transverse momentum ranging from 
20-320 GeV (corresponding to ISAJET's PT and QTW cards). The Monte Carlo rou- 
tine inefficiently samples over large j?t ranges, so each calculation is automatically 
broken into subranges of 20-40, 40-80, 80-160 and 160-320 GeV, and the results 
combined. We typically generate 50 fb -1 of sample for each signal and background 
case. 

3.2.2 Detector simulation 

As each event is generated, it is processed through a toy detector simulator consisting 
of the following elements (where ij is the pseudorapidity, cp is the azimuthal angle, 
and AR = ^/(Ar)) 2 + (A</>) 2 ). 

• Calorimeter simulator: We implement a toy calorimeter based on the IS A JET 
CALSIM subroutine. The segmentation is Ar] x A<f> = 0.1 x 0.087 extending 
to a rapidity of \rj\ = 4. There is a hadronic calorimeter, into which hadrons 
deposit their energy with a resolution of 50%/^/i? T , and an electromagnetic 
calorimeter which captures electrons and photons with resolution 15%/y / £' T . 
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We do not attempt to simulate effects of cracks or dead regions that are specific 
to particular detectors. 

• Isolated lepton identification: We sum the hadronic transverse energy in a 
cone of AR < 0.4 around each lepton. If this hadronic energy is less than 
25% of the lepton's transverse energy, then the lepton is declared isolated. p T 
thresholds for isolated leptons are given for the relevant channels. 

• Jet identification: We use ISAJET's GETJET jet-finding algorithm. Jets are 
defined as hadronic clusters with total E T > 15 GeV falling within a cone 
of radius AR < 0.7 and subject to \tj\ < 3.5. We do not perform jet energy 
corrections. 

• Silicon vertex detector (SVX): We simulate a SVX detector for tagging 6-jets. 
We identify each weakly-decaying B hadron in an event with E T > 15 GeV 
and \r]\ < 2. If AR(B, jet) < 0.5 for some jet then that jet is tagged as a 6-jet 
with an efficiency of 50%. 

• Charge multiplicity counter: The charge multiplicity of a jet is defined as 
the number of charged long-lived hadrons with E T > 0.5 GeV lying within 
AR < 0.5 of the jet axis. We use this for r veto studies. 

• Top mass reconstruction: For single-lepton events we perform tt event recon- 
struction and calculate the resulting top mass. Details of this procedure are 
given in Chapter |5[ 

Events are produced by ISAJET, passed through this detector simulator and 
subjected to software "trigger" cuts. Those which pass these cuts have their event 
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data written. This data is then post-processed to prepare ntuple files for later cut 
analysis using the CERN Physics Analysis Workstation (PAW) software^]. 

3.3 Methodology 

We study the decay modes t\ — * bW± and t\ — > cZ\ separately, in each case taking 
the t\ branching fraction to be 100% for the given mode. For the t\ — > cZ\ mode, 
the stop signal then depends only on the two SUSY masses and Tn^, and is 
independent of mixing angles in the top squark and gaugino sectors. 

On the other hand, signals from t\ — ► bW\ mode events (with the cascade 
decay W\ — ► Z\ -l-fermions) are determined by the three SUSY masses , and 
, together with the branching fractions for W\ — * e, fj,, t. In most of our analysis, 
we collapse this parameter space from four dimensions to two. First, we generally 
take m~ = /2, which is approximately true in many interesting classes of 
SUSY models (see eqn. ( |1.6| )). This only requires that the lighter chargino and the 
two lighter neutralino are gaugino-like, and that the gaugino masses obey the mass 
unification condition. While this occurs canonically in grand unified models, there 
are also many other models where this condition is obeyed but for entirely different 
reasons. Second, we set the chargino's leptonic branching fraction equal to the Ws 
leptonic branching fraction: 

B(W! -»• e) = B(Wi -»/i) = B(Wi -» r) = B(W -»• e) » 1/9, (3.2) 

as discussed earlier. We also examine how our conclusions are affected when these 
restrictions are relaxed. 

For each channel investigated, we generate event data sets for a grid over 
the appropriate parameter space of relevant SUSY masses, together with event data 
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sets for the relevant background processes. The grid spacing is usually Ami = 10 
GeV, Am^ = 10 GeV and Am- = 5 GeV. These data sets are then subjected to 
a variety of cuts on measurable quantities designed to maximize the discovery reach 
over the SUSY parameter space, using the PAW software to analyze cuts. 

We consider a particular SUSY model "discoverable" if (1) there is a 5a 
statistical significance 

vm > 5 ' (3 - 3) 

where N s (iVj,) is the expected number of signal (background) events, (2) the ex- 
pected signal level is at least 25% of background, 

^ > 25%, (3.4) 

and (3) at least 5 signal events are expected 

N s > 5. (3.5) 

When comparing our results to those of other studies, bear in mind that many au- 
thors report their results at a 3a level rather than our 5a level. The 25% requirement 
eqn. (|3.4j ) encodes our overall uncertainty in the calculated absolute cross-sections 
for background processes. And we require that at least 5 signal events be expected 
for protection against uncontrolled non-physics backgrounds. 

A signal which fails N s /Ni, > 25% is called "background limited." Increasing 
the integrated luminosity C cannot help this situation, and only clever cuts can 
rescue the signal. When a signal fails N s /^/N^ l > 5 it is "rate limited." A rate- 
limited signal that just satisfies N s /Nb = 25% can be recovered by increasing £ 



until the point N s /y/Nb = 5a is reached. At this point, where the equalities in (3.3) 
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and (0 both hold, we have N B = 400 and N s = 100. Writing N B = a B C, with 
a B the background cross-section, then 

L = when = 5a and — = 25% (3.6) 

<j b vm N b K J 

which will be useful to us in the sequel. 

The observability criteria adopted here are conservative, and elaborate sta- 
tistical schemes might well extend our reported reach considerably. We have, for 
instance, not investigated the utility of using information from the background-rich 
kinematical regions excluded by our cuts to get a better data-derived background 
estimation which would allow us to relax our 25% criterion. Also, it is worth re- 
membering that by combining results from both CDF and DO one can have get a 
stronger signal than either detector would get alone. 
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Chapter 4 

The decay mode t\ — ► cZ\ 



If the top squark decays via t\ — > c-Zi, it can be searched for in the multijet + 
$r channel, which is the standard channel for SUSY searches at hadron colliders. 
This channel picks out events where the two Z\ neutralinos escape detection in the 
experimental apparatus, resulting in a large net missing transverse energy. We also 
investigate the possibility of tagging one of the charm jets via its muon decay to 
further enhance the signal over SM backgrounds. We take the branching fraction for 
t\ — > cZ\ to be 100%. The SUSY signal in each of these channels is then completely 
determined by the two quantities and m^r . The stop mass sets the produc- 
tion cross section and the two masses together determine the kinematics. We use 
ISAJET to simulate signal events over a grid in the — plane. For each value 
of (m^ , fn~^) we generate events, run them through our toy detector simulation and 
compare the resulting cross-sections with those from the SM background processes 
to determine the reach of the Tevatron experiments. 
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c jet 




c jet 



Figure 4.1: The standard event topology for the^x + jets channel of the t\ — > cZ\ decay 
mode. 

4.1 Missing E^ + jets channel 

The t\t\ event topology is exhibited in Fig. |4.1| . The two LSPs combine to provide 
a large and the two charm quarks are required to provide hard jets. Standard 
Model events that have hard neutrinos (to supply ^-p) without isolated leptons, as 
in Fig. are the main physics backgrounds to this t\ — ► cZ\ search. Z — » vv 
is a problem when direct high p T weak boson production is accompanied by gluon 
or quark radiation, and also when an outgoing quark from a QCD hard-scattering 
radiates a Z which decays invisibly. We use jet hardness and geometrical cuts to 
suppress these relative to the signal. Similarly, events with W — » tv contribute to the 
background when the r decays hadronically or into a non-isolated lepton and fails to 
be distinguished from a hadronic jet. The processes of concern are qq' — > W — > tv 
plus one or two QCD jets and gq — > Wq' — > rvq' . Ws from ti events also contribute 
to the background. 

Note that within IS A JET, Z-production processes such as gq — > qZ with final- 
state gluon radiation (Fig. |4.3| a) are treated separately from QCD hard-scattering 
processes such as gq — ► gq or q'q — > q'q with final-state Z radiation (Fig. |4.3|b). 
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While the former are readily and efficiently simulated with the ISAJET DRELLYAN 
reaction card, the latter are treated by ISAJET as QCD processes and it is difficult 
to obtain a large Monte Carlo integrated luminosity for the latter due to the huge 
QCD hard scattering cross section. By analyzing ~ 10 7 events, we found that these 
radiative events add about 5% to the related qq — > gZ background. 

4.1.1 Run I 

We performed an early analysis to estimate discovery prospects in this channel 
with the Tevatron Run I data. Based on this analysis we recommended a set of cuts: 
(1) $ T > 50 GeV to reduce backgrounds from QCD heavy flavors and mismeasured 
jets; (2) at least two jets in the event, one of which is central (\r/\ < 1); (3) all jets 
separated by at least 30° in azimuth from^ T , and a jet-jet separation less than 150° 
in case there are only two jets; and (3) no identified leptons. We also prescribed a 
compound cut (4) p T (fast jet) > 80 GeV for A(f)(]fi T ,j) > 90°, else p T (fast jet) > 50 
GeV. We found the dominant background after these cuts to be W — > tv, with the 
r decaying hadronically. Depending on the ability of experimentalists to reject the 
W — ► tv background and subtract the Z — » vi> backgrounds, we estimated that, 



t jet quark jet 




gluon jet ISR jet 



Figure 4.2: Some representative background processes for the-^T + jets channel. 
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Figure 4.3: Examples of ISAJET Z radiative processes. Diagram (a) is a well-controlled 
ISAJET DRELLYAN process; processes like (b) require simulation of a very large number of 
ISAJET QCD events. 

depending on m-x , stops as heavy as 100-125 GeV might be probed in this channel 
with 100 pb _1 of data. 

So far, the full Run I data set has not been analyzed in this channel. The 
DO Collaboration f§ looked at 13.5 pb" 1 of data from the 1992-1993 run. They 
required $ T > 40 GeV, two jets with p T > 30 GeV, 90° < A(j)(jx,j 2 ) < 165°, 10° < 
A0(^ T ,j), A0(|? T , j) < 125° and no identified leptons with p T > 10 GeV. They 
excluded a region at the 95% confidence level of ^ 90 GeV, — ^ 35 
GeV. 

4.1.2 Run II 

The larger data sample anticipated in Run II enables us to sharpen the selection 
cuts for the stop signal. For the^ T + jets channel we require the following initial 
cuts: 

(i) $ T > 50 GeV; 

(ii) at least two jets with p T > 50 GeV, one with \rj\ < 1; 
(Hi) for all jets, A</>$ T ,j) > 30°, 

and if only two jets then A(ft(ji, j 2 ) < 150°; 
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(iv) no identified isolated e or //; 
(v) no SVX b tag. 

These cuts will reject QCD backgrounds from heavy flavor production (g — > cc, bb) 
and multijet production with mismeasured jets. Cuts (i) or (ii) may serve as trig- 
gers. Note that we require very hard jets, one of which must be central. Cut (Hi) 
is designed to reduce detector-dependent backgrounds due to mismeasured QCD 
jets. Cut (iv) suppresses background events with leptons from W — > Iv, £ = e, /i. 
Cut (v) suppresses the ti background by vetoing events with 6-jets which will be 
efficiently tagged by the Run II Tevatron detectors. (In our simulation we ignore c 
contamination of SVX tags. A closer analysis which pays attention to the detailed 
experimental properties of the particular vertex detector, including a scheme for the 
outright identification of displaced charm vertices, could improve the signal yield in 
this channel. See the remarks at the end of the next section.) The cross sections 
for the backgrounds and selected SUSY signal cases after these cuts are displayed 



in Table iA in the column labelled (v). 

For the Z — > vv background due to Drell-Yan qq annihilation, the two hard 
jets come from QCD radiation recoiling against the Z boson; these jets are often 
close together in azimuth. The c-jets in the signal, on the other hand, come from 
the two t±s which recoil mainly against each other, and hence tend to be more back 
to back in (p. Events from the W — > tv background exhibit a spectrum intermediate 
between these two; when the two hardest jets come from QCD radiation then they 
are generally close, as with Z — > vi>, but when one of the hard jets is from a 
hadronically decaying r, then the jets can have a large separation. Distributions of 



the azimuthal separation of the two hardest jets, A0(ji, j'2), are displayed in Fig. 4.4 
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background 
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Z — ► vv 


484 
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24.4 


ti (175) 


32 


24 


24 


7 


2.5 


total 


1016 
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40.3 


25% 
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92 


73 
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(vi) 
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146 
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47 


12.1 


150, 80 


146 
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86 


58 


19.7 


180, 100 


68 


54 


42 
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Table 4.1: $r + jets channel cross-sections, in fb, after cuts as discussed in the text. 
The line labelled "25%" shows 25% of the total background, and the "5<r" line shows 
the signal level needed to produce a 5a excess at 2 fb~ l . (The parenthesized 5a 
value in the last column is for an integrated luminosity of 25 fb .) For the signal 
cases, t\ and Z\ masses are given in GeV. 

for the W and Z backgrounds and a signal case of a 150 GeV t\ decaying to an 80 
GeV Z\. Imposing the cut 

(vi) A0(ji,j 2 )>9O 

eliminates about half of the Z — ► vv background along with 20% of the signal and 
other backgrounds. 



We see from Table 4.1 that while a 5a signal is already possible even for 



mi 



150 GeV, the signal to background ratio is less than 25%. To enhance the 



signal further, we note that these signal events generally have harder jets than the 
W and Z events, except when the t\ — Z\ mass difference becomes small. Also, 
the t\t\ events may be accompanied by QCD jets. The distribution of total jet 
transverse energy 



J T = ]T£ T (jet) 

jets 



(4.1) 
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Figure 4.4: Distribution of A<j>(ji,j2) in the $t + jets channel for a signal case of 150 
GcV t\ — ► 80 GcV Zi and backgrounds W tv and Z — » z/z/. The vertical axis on the 
main plot is in fb/bin. The top line in the histogram plot shows the sum of signal plus 
both backgrounds. The small subplot at the top shows the ratio of the normalized signal 
distribution to the normalized signal+background distribution. The box heights in this 
subplot are proportional to their calculated uncertainties. 
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Figure 4.5: Distribution of Jt in the^T + jets channel for 150 GeV t\ — > 80 GeV Z\ and 



backgrounds W^tv and Z—n>v, after applying cut (vi). Figure elements are as in Fig. 4.4 



exhibited in Fig. 4.5 reflects this. We make the cut 

(vii) J T > 175 GeV, 

which reduces the W and Z backgrounds by 40% at a cost of a fifth to a third of 
the signal. 

Further cuts based on simple kinematical quantities do not seem to be avail- 
able. In order to extend the useful discovery reach in this channel we investigate 
the possibility of discriminating the rs in the W — > rv background. Since we have 
already vetoed events in this channel with an identified lepton, we expect most of 
these rs to decay hadronically, producing jets with low charge multiplicity. For the 
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purposes of our simulation, we calculate the charge multiplicity of a jet, mult(j), 
as the number of long-lived charged hadrons with Ej, > 0.5 GeV lying in the jet 
cone. Distributions of the minimum charge multiplicity, mm,- mult(j), are shown in 
Fig. [4.6[ Bare ts will canonically decay to one or three charged hadrons, but upon 
analyzing simulated events we find that a "r-veto" requirement 



mhij mult(j) > 4 



gives the best separation of signal and background over the SUSY parameter space, 
eliminating almost 60% of the background along with 30-40% of the signal. Note 
that 70% of the tt background also fails this cut since t decays with large ^ T (from 
an energetic u) and no visible leptons will usually involve W — > tv. 

The results of applying cuts (i)-(tim) are displayed as a contour plot in 



Fig. 4.7. This plot shows the SUSY parameter space m^-m^ . The two diagonal 
lines indicate the kinematical boundaries of the decay t\ — > cZ\\ the upper left hand 
corner is excluded by the requirement m~ < that Z\ be the LSP, and in the 
lower right corner the 3-body decay t\ — > bW Z\ opens up. The combined LEP 
II/Tevatron Run I excluded area (see Fig. |2.6|) is denoted by the hatched region. 
The legend in the upper left corner shows the estimated background cross-sections, 
in fb, after cut (viii). The contour labels are also given in fb. The heavy solid line at 
48 fb represents N s /^/T7h = 5a for the Run II integrated luminosity of 2 fb -1 . The 
heavy dashed line is at 45 fb = 25% of the 180 fb total background. From the plot 
we see that even with our conservative criteria, experiments at the Tevatron can 
discover most cases with ^ 70 — 80 GeV and ?S 165 GeV in this channel. 



Note that the contours of Fig. 4.7 extend into the region marked t\ — > bWZi, 



where the three-body decay is kinematically open. Above the t\ — > bW Z\ line the 
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Figure 4.6: Distribution of minj mult(j) in the + jets channel for 150 GeV t\ — > 80 
GeV Zi and backgrounds W — > tv and Z — > vv, after applying cut (vii). Figure elements 



are as in Fig. 4.4 
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m( t,) (GeV) 

Figure 4.7: Cross-section contours, in fb, for the fl T + jets channel after cut (viii) as 
described in the text. The heavy solid line at 48 fb is the 5<r discovery limit for an integrated 
luminosity of 2 fb -1 . The dashed line is 45 fb of signal, which is 25% of background. The dot- 
dashed line shows the reach at 25 fb -1 after cut (ix). The hatched region is experimentally 
excluded at the 95% confidence level. 
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branching fraction B(t\ — > cZ\) is taken to be 100%, while below the line the two 



decay modes are expected to compete |87|, p6[ . 
4.1.3 Run 11+ 

This channel is background limited in that even for higher integrated luminosity 
the search is limited by the signal to background ratio, and without refining the 
requirements on event topology, such as the c-tagging studied in the next section, 
little improvement seems available for the high luminosity scenario of the Run 11+ 
Tevatron. The availability of more events, though, allows us to make deeper cuts 
and still preserve an acceptable N s /y/Nb. By tightening the previous cuts to 

K) > 120 , 

(vii') J T > 225 GeV, 
[viii') minj mult(j) > 5 

we can push the discovery limits at the heavy t\ end of the parameter space, up to 
about 185 GeV. Our projection of the reach at 25 fb _1 after these cuts is shown 



on the contour plot of Fig. 4.7 as a heavy dot-dashed line. We see that the reach 



extends beyond mx = 180 GeV. 

4.2 c-tagged channel 

In the c-tagged channel we attempt to identify the charm quark by looking for a jet 
with a soft muon due to the decay c — * sjxv^. The backgrounds here are similar to 
those in the previous channel. Z —* vv contributes when accompanied by tagged 
heavy flavor production (e.g., g — ► cc from initial state QCD radiation). W — ► e,fj, 
and W — > t — > e, /i are a problem when the lepton mis-tags an unrelated jet. ti 
events will be a reducible background when a semileptonic 6-jet is tagged as a c-jet. 
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4.2.1 Run I 

Our earlier analysis of this channel geared for the Run I |2q| integrated luminosity 
suggests cuts as in the Run I$ T + jets channel (see Sect. |4.1.1|) , except that we here 
require a tagging muon with p T (/^) > 3 GeV within AR = 0.4 of a jet axis, and the 
compound cut (4) is changed to (4') A0($? T , nearest jet) < 90° or p T (fast jet) > 50 
GeV. We concluded that with the integrated luminosity available in Run I, the c- 
tagged channel would not extend the reach beyond the^ T + jets channel because, 
although it has a better signal to background ratio, the signal is severely rate limited. 
The CDF Collaboration has analyzed 88 pb _1 of Run I data|l9| looking for 
c-tagged jet with a different technique than ours, using the SVX detector to 
identify charmed hadrons by looking for displaced vertices. They required -$? T > 40 
GeV, 2 or 3 (hard) jets with > 15 GeV within \r]\ < 2 (to reject tt events, which 
usually have 4 or more hard jets) but no (soft) jets with E T < 15 GeV (to suppress 
QCD multijet events), 45° < A<f)($ T ,j) < 165°, 45° < A(p{j 1 ,j 2 ) < 165°, and no 
leptons with p T > 10 GeV. Then they look for the SVX charm tag. By this means, 
they get a 95% confidence level excluded region bounded roughly by mi < 115 GeV 
and m~ r < 50 GeV. 

4.2.2 Run II 

Initial cuts for the c-tagged channel are: 
{%) E T > 50 GeV; 

(ii) two jets with p T > 30 GeV, one with \r]\ < 1; 

(m) for all jets, A0(^ T , j) > 30°, 

and if only two jets then A<p(ji, j 2 ) < 150°; 
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(iv) one muon with p T > 2 GeV, \rj\ < 1.7 within AR = 0.4 of a jet; 

(w) no other visible leptons; 

(vi) no SVX b tag. 

Cuts (z) to (Hi) are implemented for the same reasons as in the previous channel. 
The charm tag lets us relax our jet p T requirement. Cut (iv) identifies the muon 
tagging a c-jet. After this cut, some background events will have genuine c-jets, 
others will have /Lt-tagged b quarks, and yet others will get fake c-tags due to mis- 
cellaneous muons accidentally coinciding with unrelated jets. Most ti events in this 
sample are of the b — > fi sort. Significant numbers of Z — » vv events are also of this 
type, since Z — > vv events which pass the initial cuts are almost all due to gluon 
splitting to heavy flavors, and g — > bb occurs roughly as often as g — > cc for these 
events after the hard cuts (i) and (ii). We veto extra isolated leptons and identified 
6-jets with cuts (v) and (vi). 

After these cuts, we still have the reducible background from non-SVX-tagged 
6-jet production as well as fake c-tags, in addition to the irreducible background with 
real c-tags. To reduce the former, we consider the transverse momentum of the muon 
with respect to its tagged jet. This is 

Pt 1 = \M x PU)\ » 

where p(j) = p(j)/\p(j)\ is the jet direction unit vector. Since m& > m c , the average 
pip / wm be greater for 6-jet samples than for c-jet samples. Figure [4.8| displays 
scatter plots of versus the /U-jet axis separation AR(fi,j) for the backgrounds 
together with a pair of representative signal samples. Also shown is the line AR + 
pf/h GeV = 0.4 which we have selected for our cut. The efficient separation is 
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Figure 4.8: Scatter plots of p^ 1 vs AR(/i,j) for backgrounds and selected signal points in 
the c-tagged channel after cut (vi). The diagonal line shows the cut (vii) line AR + p^ l /5 
GcV=0.4. 
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Table 4.2: c-tagged channel cross-sections, in fb, after cuts as discussed in the text. 
The line labelled "25%" shows 25% of the total background, and the "5cr" line shows 
the signal level needed to produce a 5a excess at 2 fb^ 1 . The parenthesized 5a value 
in the last column is for an integrated luminosity of 25 fb . For the signal cases, 
t\ and Z\ masses are given in GeV 

obvious from the diagram, and we enhance the c-jet purity by requiring 
(vii) AR(nJ) + p T T el (fi\j)/5 GeV < 0.4. 



As shown in Table 4.2, this cut eliminates 60% of the background at a typical cost 
of 20% in signal. The large reduction shown in the table for W — * Iv is due to 
the limit on p 1 ^, since the high-p T muons in these events are only coincidentally 
associated with a jet. We have checked further cuts, such as limits on p T (£), and 
found that they have nothing substantial to add after cut (vii). 

A contour plot summarizing the result of applying cuts (i)-(vii) is shown 



in Fig. 4.£. The total background remaining at this level is 52.2 fb. The heavy 
solid line in the figure gives the 25.5 fb signal level, which corresponds to 5a at an 
integrated luminosity of 2 fb -1 . This channel affords substantially more coverage 
than the^ T + jets channel of the previous section, and much of the parameter space 
m; ^ 170 GeV, m-x ^ 90 GeV is accessible to detection. 
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Figure 4.9: Cross-section contours, in fb, for the c-tagged channel after cut (vii) as 
described in the text. The heavy solid line at 25.5 fb is the 5cr discovery limit for an inte- 
grated luminosity of 2 fb -1 . The dashed line is 13 fb of signal, which is 25% of background. 
The dot-dashed line shows the reach at 25 fb -1 after cut (viii). The hatched region is 
experimentally excluded at the 95% confidence level. 
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Some recent analysesjlS, 51] have looked at tagging charm jets directly with 
the SVX. The CDF Collaboration has done this for the Run I data, yielding the 



excluded regions shown on Figs. |4.7| and |4.9| , They rely on improved resolution to 
identify the small impact parameters of c-decay displaced vertices. This analysis has 



been extended for the Run II scenario [51]. These searches proceed in a complemen- 
tary channel to the one we investigate here, since they accept SVX tags and reject 
identified leptons (to reduce W and Z backgrounds) while we veto SVX tags and 
require a lepton. The CDF estimated reach with these cuts is similar in extent to 
ours. It may be worth investigating whether combining the two c-tagging techniques 
would further extend the Tevatron reach. 



4.2.3 Run 11+ 

Unlike the^ T + jets channel, the signal in the c-tagged channel is not background 
limited at 2 fb -1 , as indicated by the large gap between the solid 5cj line at 25.5 fb 
and the dashed N s /Ni, = 25% line at 13 fb. Thus, more integrated luminosity can 
rapidly improve the reach with these same cuts. By eqn. fl3.6| ), 7.7 fb -1 will bring 
us to the 13 fb contour. After this point further cuts are needed for the Run 11+ 
scenario. 

To cut the background further, we note that by vetoing events with non- 
tagging leptons we have increased the fraction of hadronically decaying is, so that 
most tt events surviving cut (vii) are either t —* W — > r decays or Iv + 4-jet 
events with the lepton faking a charm tag. About 60% of these events have 5 or 
more reconstructed jets, as seen in Fig. 4.10) , so limiting the jet multiplicity helps 



suppress tt. 
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Figure 4.10: Distributions of A<P(c,$t) and nj e t s in the c-tagged channel for a signal case 
of 180 GeV h -> 110 GeV Z 1 and backgrounds it, Z ->vD,W -> tv and W -> tu, after cut 
(uii). The vertical axes are in fb/bin. 
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As in the previous channel, the jets in Z — > vv recoil against the Z boson, 
which itself is converted to^ T . Thus the angle A(/>(c,-^ T ) between the tagged jet 
and Jfj^ will be large. For t\t\ there are two c-jets recoiling independently against 
their two LSPs, and the distribution of Ac^(c,^ T ) should be broader. This is indeed 



the case, as shown in Fig. 4.1C. In W — > Iv events with p T (W) <C mw, the angle 
A<p (I, u) will be large and if I mistags a jet and v supplies ^ T then A0(c,^ T ) will 
be large too. We address these points by insisting that $ T and the c-jet not be too 
back-to-back (75% of the Drell-Yan background has A0(c,^ T ) greater than 120°). 
Our Run 11+ optimized high luminosity cut is 

(viii) A(f)(c,$ T )< 145° and n jets < 5 

which eliminates over half of the background events and fewer than 1/3 of the 



signal events. Cross-sections after this cut are shown in Table 4.2. The dot-dashed 



discovery line in Fig. 4.9 indicates that t\s as heavy as 200 GeV may be accessible 
to the Run 11+ Tevatron in this channel, even if the Z\ mass is well beyond 100 
GeV. 



It is amusing to note, with reference to Fig. O, that if (1) the electroweak 
baryogenesis story of Section [2.1.1| is correct and mi < 160 GeV, and (2) the galactic 
halo story of Section |2.1.1| is correct and m^ — m~ > 20 GeV, and (3) we "double" 
the integrated luminosity by combining the CDF and DO experimental results, and 
(4) the stop decays as t\ — > cZ\ without the 4-body decay of Section [2]^ having a 
significant branching fraction, then the stop should be discovered at the Run 11+ 
Tevatron! 
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Chapter 5 

The decay mode t\ — > bW\ 



When mi x > nib + m ^ 1 i the stop decays via t% — > bW\ with essentially 100% 
branching fraction since competing decays occur only at higher order. The W\ 
then cascades to a Z\ plus a SM fermion pair. Unlike the t\ — > cZ\ mode, where 
the analysis depended only on and m~ , here we have 4 parameters: the three 
masses , , m~ and the leptonic branching fraction of the chargino, B(W± — ► 
e, fi or t) (assuming lepton universality in W± decays). To make the parameter space 
tractable, we will adopt the gaugino unification relation w Imj. (eqn. (|1.6|) ) 
for most of the discussion, and separately examine the consequences of relaxing this. 
We also will take B(Wi — ► ~ ^(VF — > £) for the main analysis. 

When both WVs decay hadronically, the signal is jets plus^ T from the pair 
of Zis. This signal, with degraded ^ T due to the two-step decay of the tx, is difficult 
to discriminate from SM processes. Therefore, we focus on channels where one or 
both Wi's decay leptonically. In the case of one leptonic decay we look for a hard 
isolated lepton together with a 6-tag along with the usual^ T . We refer to this search 
as the "6-jet + lepton channel". In the "dilepton channel" we look for ^ T plus 2 
unlike-sign isolated leptons (e or u) as the result of both W\s decaying to leptons. 

60 



The dominant background in both channels is SM ti production and decay 
followed by the leptonic decay of one or both W's; such processes have 6-jets, leptons 
and substantial $ T from the associated neutrinos. Indeed, the decays t —> bW — > 
bff and t\ — ► bW\ — > bj fZ\ differ only in the presence of two Z\s in i\i\ events 
and the presence of on-shell Ws in ti events; event topologies from stop and top 
production are identical. The two massive ZiS in the final state lead to a distinctive 
softening of the spectra of many kinematical quantities in the signal events. The 
on-shell W in t — ► bW events implies useful kinematical constraints on its daughters 
/,/'. This is especially helpful when only one of the Ws decays leptonically, in 
which case the transverse mass of the lepton and $ T help to discern the real W in 
the event. The W kinematics are also useful in suppressing other backgrounds, such 
as those due to Drell-Yan Ws, W radiation, and W pair production. 

5.1 6-jet + lepton channel 

In the 6-jet + lepton channel, we look for^ T , hard jets with a tagged b, and an 
isolated e or fi. In this channel we contend with two main SM backgrounds. The first 
is ti production, as mentioned above. The second is due to events with a leptonically 
decaying weak boson W — > e, fi, r associated with QCD radiation. This reaction has 
a huge cross-section, and readily delivers large ^ T from the decay neutrino, along 
with a hard isolated lepton. Much of this background is suppressed by requiring a 
tagged 6-jet in the event. The W jets that pass this criterion generally have 6s from 
gluon splitting g — > 66. A minor background comes from the process Z — > r + r~ 
when one of the rs decays leptonically and a real or fake 6 is tagged. We ignore 
mistagged 6-jet backgrounds. 
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5.1.1 Run I 



In our earlier pre-LEP analysis [28] of the Run I situation for this channel, we 
suggested the following initial cuts: (1) $ T > 25 GeV; (2) at least 2 jets with 
p T > 15 GeV, one of which lies in the central region \rj\ < 2; (3) an isolated electron 
or muon with p T (e) > 10 GeV and p T (fi) > 5 GeV; and (4) an SVX 6-tag (with 
an efficiency estimated at 40% (in the barrel) x 30% (identified displaced vertex) 
= 13%)Q The purpose of these cuts is the same as the analogous cuts in our Run 
II analysis presented below. After these, we also recommended (5) no more than 
4 jets (with p T > 15 GeV); and (6) the transverse mass cut m T (l,]fi T ) < 45 GeV 
(see below for this quantity). We found that in this channel the signal dwindles as 
mjj — ► rriyy , where the b becomes too soft to tag. We determined that the region 
Tifj ~ 90 GeV, tu^ ;$ 60 GeV would be accessible to this treatment with the Run 
I integrated luminosity of 100 pb _1 . 

The CDF Collaboration examined this channel for 88 db 4 pb _1 of their Run 
I data@. They used the cuts (1) $ T > 25 GeV; (2) at least 2 jets with p T > 12 
GeV for the hardest jet and p T > 8 GeV for the second jet; (3) an isolated electron 
or muon with p T > 10 GeV; and an SVX 6-tag. They further required that (5) 
any dilepton have invariant mass m (£ + £~) < 60 GeV; and (6) the azimuthal angles 
between ^ T and each of the two hardest jets satisfy Ac/> (|£ T , jet) > 0.5 rad. They 
found no evidence for t\t\ production, but were not able to improve the existing 
(LEP II) stop mass limits. 



x In the meantime, SVX technology has improved considerably. Resolution has 
gone up and coverage has extended to \tj\ < 2. For Run II, both DO and CDF will 
be equipped with SVX detectors. 
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5.1.2 Run II 

For the higher luminosity run, our initial cuts for this channel are 
00 $ T > 25 GeV; 

(ii) exactly one isolated e or fi with p T > 8 GeV and \t]\ < 2: 

(Hi) at least two jets with p T > 15 GeV and |r/| < 2; 

(iv) at least one SVX 6-tagged jet with \rj B \ < 1.4. 

Cuts (i) to (Hi) are canonical for this channel. The 6-tag, cut (iv), is critical to 
reduce the W background. The 6-tagging efficiency for signal events ranges from 
one third to two thirds, while less than a percent of the W events get tagged. The 
\rj B \ distribution for W events with 6s due to g — > bb is broader than that for signal 



events, as shown in Fig. 5.1, hence the cutoff at \t] b \ < 1.4. The soft-lepton 6-tagging 
used in SM ti studies is avoided here to avoid the large background from events like 
sg — > Wc — » ivc or qq — > Wg — > tvdcc where t is isolated, v provides ^ T , the c is 
tagged and g provides the second jet. (In the ti studies, one makes hard cuts to 
reject these W events; as we shall see below such hard cuts would remove too much 
of our stop signal.) 

We expect that for ti and W — > e, [i backgrounds, where missing transverse 
energy is due to a single neutrino and^ T ~ p T (y), the transverse mass m T (^,^ T ) 
of the lepton and^ T should show a strong Jacobian peak near my/, while the signal 



distribution should be a broad bump. Figure 5.2 shows that this is indeed the case. 
The transverse mass is the invariant mass of the^ T vector and the projection of the 
lepton's momentum in the transverse plane, 

(l$ T f = 2p T (£)$ T (1 - cos Acp) . 
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Figure 5.1: Distribution of \t]b\ in the 6-jet + lepton channel for the initial cut sample. 
The signal case is a 160 GeV ti decaying to a 120 GeV W\ which in turn decays to a 60 
GeV Z\. The vertical axes are in fb/bin. 
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Table 5.1: 6-jet + lepton channel cross-sections, in fb, after cuts as discussed in the 
text. Signal levels required for detection are given in the "25%" and "5<r" rows. 
The parenthesized 5a value in the last column is for an integrated luminosity of 25 
fb -1 . For the signal cases, t\ and W\ masses are given in GeV, and m~ = /2. 
Cut (ix) is made after cut (viW). 

We can get even more information from this kinematical data by forming a "semi- 
transverse" mass which is the invariant mass of-^ T and the lepton's full 4-momentum: 



m T+ (^,$ T ) 2 = 2p T (t)$i (1/ sva.'di — cosA^j) 



(5.1) 



where •&£ = 2 arctanexp (— rji) is the lepton's polar angle. The m T+ distribution is 
also plotted in Fig. |5.2| , where one sees that for the VF-containing backgrounds it 
has a broader high-end tail than m T , and that the distribution vanishes in the limit 
m T+ — > 0, unlike m T which remains finite there. Both of these behaviors make m T+ 
a better discriminator than 771/ rp j QtS the signal distributions for the two quantities 
are largely the same. With a cut of 

(v) m T+ (£,# T )< 60 GeV 

3/4 of the tt and W — > iv background is removed with less than half of the signal 
lost. (In comparison, an m T cut removes less than 2/3 of the background Ws at the 



same signal loss.) The results of this cut are listed in Table 5.1 
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Figure 5.2: Distributions of tot(^,^t) and mT+(f,|T) in the b-jet + lepton channel, after 
cut (if). The signal case is a 160 GeV t\ decaying to a 120 GeV W\ which in turn decays 
to a 60 GeV Z\. The vertical axes are in fb/bin. The dashed lines show cut (v). 
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The ti background still dominates the signal. To reduce it further, we note 
that the lepton from W decay is typically more energetic than that from W\ decay. 
In the decaying W or W\ frame, one can crudely compare the available lepton energy, 
mpi//2 for W — > £v, with (wt^ — m~^)/2 ~ ra wJ^ L ^ or ^ ~~ * anc ^ use ^ ne ^ ac * 
that generally < 2miy for the models accessible in this channel. We have also 
found that the third jet tends to be harder in ti events than in our signal events. 
The reason is similar to that given for the lepton since, in the absence of hard QCD 
radiation, this third jet mainly comes from the hadronically decaying W or W\. We 
capture this feature by defining the quantity 



W = ^fp T (l) 2 +P T (h) 2 (5.2) 
where p^(js) is taken as for those events with only two jets. Imposing the cut 

(vi) W < 60 GeV 

reduces the dominant background from tt by almost half and the signal by just 



10-20%. W distributions are displayed in Fig. 5.3, and the results of the cut are 



tabulated in Table 5.1 



We turn our attention now to the W events, which are distinguished in several 
ways from the signal (and tt) events. First, their jets come from QCD radiation so 
they have a lower jet multiplicity. Also, they have a low efficiency for a second b tag, 
a fact used in the Run 11+ section below. Finally, the W and Z events are generally 
softer than the signal events, which is reflected in the distribution of the quantity 

H T+ =P T+PT (e) + J2 PT (j), (5.3) 
j 

which is well-known from ti studies to provide a good separation of ti from the 



W — ► e, /i, t background. Figure 5.4 shows these distributions. The t\ distribution 
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Figure 5.3: Distribution of W in the 6-jet + lepton channel, after cut (v). The signal case 
is a 160 GeV t\ decaying to a 120 GeV W\ which in turn decays to a 60 GeV Z\. The 
vertical axes are in fb/bin. The dashed lines at W = 60 GeV indicate cut (vi). 
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unfortunately falls in between the two backgrounds, a predicament which is quite 
general across the parameter space we have investigated for this channel. Still, the 
cut 



(vii) H T+ >140 GeV 
eliminates 40% of the W background with only a few percent loss of signal, as shown 



in Table 5.1 



So far, we have not used the fact that elaborate techniques have been devel- 
oped for reconstructing tt events. The single lepton plus 6-tag channel is a primary 
source of tt events for Standard Model top mass reconstruction studies and our 
experimental colleagues will certainly subject events in this channel to a thorough 
characterization using these techniques. Mass reconstructions begin by assuming 
that the missing Er represents the momentum of the v from a leptonic W — > Iv 
decay (the other W is assumed to have decayed hadronically) . By forcing v and t to 
reconstruct a real W, one gets two solutions for the longitudinal momentum of the 
v. Together with combinatorics from assigning jets to the (possibly tagged) bs, the 
hadronic W daughters and miscellaneous QCD radiation, one gets several fitted so- 
lutions. For top mass reconstruction studies, these solutions are variously all used, 
used with weights, or judged to select a best solution. Although straightforward 
in principle, such analyses are quite complicated in practice. Results are sensitive 
functions of jet energy corrections, and thus depend strongly on the details of the 
detector simulation. 

To estimate the utility of using tt characterization to suppress these events as 
a background, we have performed toy reconstructions on our Monte Carlo samples. 
We use the procedure outlined above, without applying jet energy corrections. For 
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Figure 5.4: Distribution of Ht+ in the 6-jet + lepton channel after cut (vi). The signal 
case is a 160 GeV t\ decaying to a 120 GeV W\ which in turn decays to a 60 GeV Z\. The 
vertical axes are in fb/bin. The dashed line at ifr+ = 140 GeV shows cut (vii). 
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Figure 5.5: Distribution of fh t in the 6-jet + lepton channel after cut (vii). The signal 
case is a 160 GeV i\ decaying to a 120 GeV W\ which in turn decays to a 60 GeV Z\. The 
vertical axes are in fb/bin. Events that are not reconstructible (for which fht = 0) are not 
shown in the plots; the plotted events constitute 82% of the tt sample, 15% of the W ' — > Iv 
background, and 46% of the signal sample. 
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each event, we judge a best solution to be that which minimizes the quantity 

(m lep _ m had )2 + {mw(jj) _ (5 4) 

where mj ep is the reconstructed mass of the semi-leptonic top, m had that of the 
hadronic top, and m w (jj) the invariant mass of the two jets assigned to the hadronic 
W. From this best solution we form the quantity 

fh t = (m\ ep + m^ d )/2 (5.5) 

which is our (crude) representation of the fitted top mass. Reconstruction is not 
possible for events with less than four jets, and for those events we define fht = 0. 
Events with more than four jets are reconstructed with all possible jet assignments 



to find the best fit. The resulting fht distributions are shown in Fig. 5.5. The crudity 
of our reconstruction is manifested in the tt plot, which is asymmetric and peaks 
noticeably below nit, largely due to our lack of jet energy correction. Even at this 
level of sophistication though we can reject almost half of the contamination with 
a penalty of less than 15% of signal by cutting 

(viii) fht < 150GeV. 



Figure 5J3 summarizes the discovery potential of applying cuts (i)— (viii) in 
the 6-jet + lepton channel. The parameter space of this plot is m^-m^ , where 
we have taken = and B(W\ — ► I) = B{W — ► I). The contours show 

signal cross-sections in fb, on a total background of 128 fb. The solid line at 40 fb 
indicates the 5<r reach for 2 fb -1 of integrated luminosity, which is an expected 80 
SUSY events on top of 256 SM events. The dashed line corresponds to a quarter 
of the background cross section. Stops with mass up to ~ 160 GeV should be 
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Figure 5.6: Cross-section contours, in fb, for the 6-jet + lepton channel after cut (viii) as 
described in the text. The heavy solid line at 40 fb is the 5a discovery limit for an integrated 
luminosity of 2 fb - . The dashed line is 32 fb of signal, which is 25% of background. The 
dot-dashed line shows the reach at 25 fb -1 after cut (ix). We have set m~ — for 
the signal cases. 
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discoverable in this channel for charginos lighter than 120 GeV. Notice from Table 



5.1 that for lighter stops cuts (i) - (vi) may give a better significance of the signal. 



Relaxing the assumption m~ = m^> /2 

If we consider models more general than mSUGRA and its cogeners, which have 
strict gaugino unification at the GUT scale, then the relationship = of 
eqn. |1.6| may fail to holcQ. We study this situation by performing our analysis in the 
m~ — m~ plane for the case of a 160 GeV stop. Our results are presented in the 



contour map of Fig. 5.7. The contours are signal cross-sections in fb after cut (viii) 



and we show the 32 fb and 40 fb reach lines, as in Fig. 5.6. The straight dot-dashed 



line traces = /2, and indicates the slice of this plot corresponding to the 



mi = 160 GeV slice of Fig. 5.6. The diagonal line cutting off the lower right corner 



is where Z\ fails to be the LSP, contrary to our hypothesis. 

The straight solid line marked m~ = m-x +M\y indicates where the mediat- 



ing W* in the chargino decay (see Fig. |2.5| a and attendant discussion) goes on-shell. 
The signal efficiency falls in this region as our ^-rejection cuts take hold. Also, 
when Z\ is too light, the fermions from the W\ decay are energetic enough that the 
cuts (v) and (vi) designed to suppress tt background erode the signal as well. In the 
main part of the plot, the efficiency is fairly independent of the Z\ mass, but it falls 
again in the limit m-^ — ► , as the lepton becomes undetectable. In general, the 

ti detectability declines outside the region m^/2 < < 2m^/3. 

2 Indeed, this relation is only approximately true even in the mSUGRA framework 
because of gaugino-higgsino mixing. 
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Figure 5.7: Cross-section contours, in fb, for the 6- jet + lepton channel in the m^-m-^ 
plane after cu t (v iii) as described in the text. The top squark mass is 160 GeV. Annotations 
are as in Fig. p.q. 
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Relaxing the assumption B(W± — ► £) ~ £>(VF — ► £) 

We can also ask what happens if the chargino's leptonic branching fraction is not 



VF-like. As discussed in Section [2^, this is due to /""-mediated decays, which for 
models with GUT scale universal scalar masses will usually favor slepton mediation 
(over squark mediation) due to general features of the RGE. This will increase lepton 
production and for moderate values of the enhancement 

r _ g(gW) f56) 
r " B (W -> i) (5 ' 6) 

the signal cross-sections should go up by the factor r (without any increase in back- 
grounds)]^. Of course, in the limit that W± decays only to leptons the signal will be 
lost, since we veto the second lepton in this channel. The single- lepton fraction is 
maximized when B{W\ — > £) = 0.5, or r ~ 2. If r were to actually take this value, 
then the discovery reach at 2 fb -1 would extend to the 20 fb line of Fig. 5.6, finding 
stops as heavy as 185 GeV and charginos as heavy as 140 GeV. 

5.1.3 Run 11+ 



The N S /N B = 25% line at 32 fb in Fig. |J is well beyond the N s /VNb = 5a line 



at 40 fb, which indicates that we can do more with cuts (i) — (viii) if more integrated 
luminosity is available. The 5a reach would extend out to the 32 fb line with 3.2 
fb -1 of data (see eqn. (|3.6| )), revealing stops almost as heavy as the top quark. For 
even higher integrated luminosities, we can reduce the tt background by making a 
deeper cut on the reconstructed top mass 

(viii') m t < 125 GeV. 



3 An admixture of the sfermion-mediated diagram could also change the produc- 
tion kinematics. For instance, if ~ mp t ^ m^, then the dominant ^-mediated 
decay would produce very soft leptons. 
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Then we can attack the W background by insisting on a second 6-tag. In general, 
the W event's tagged b jets are less distinguishable since they are due to gluon 
splitting (g — ► bb frequently forms only one jet). This leads to a low efficiency for 
double-tagging; while the signal and ti events have about a 25% probability of hav- 
ing a second b tag, only about 8% of the W and Z events get this second tag. We 
require a second tag in the Run 11+ scenario: 

(ix) at least 2 6-tags 

At 25 fb _1 of data these cuts give the dot-dashed line in Fig. 5.6, where the reach in 



now extends well beyond m t . The reach still cuts off at ~ 125 GeV 

though, presumably because the 6-jets have to be hard enough in order to be 
tagged. 

5.2 Dilepton channel 

In the dilepton channel we look for two opposite sign hard isolated leptons, pre- 
sumed to come from the leptonic decays of both W\S in the t\ti event. The major 
backgrounds here are ti, as usual, along with W pair production and the processes 
Z — ► £ + £~ (£ = e,/i) and Z — ► t + t~. (Note that Z — > t + t~ was left out of the 



study [28].) As in the channels previously investigated here, the ti background is 
controlled by rejecting events which are too "hard." The W pair process has a mod- 
est production cross-section, but in the absence of 6-tagging its event topology is 
quite similar to that of and this background proves the most difficult to remove. 
Drell-Yan Z — > £ + £~ has a large production cross-section, but there is no source of 
^ T in the basic reaction and those events which remain are effectively controlled by 
reconstructing the Z boson from the same-flavor dilepton invariant mass. Finally, 
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to reduce the background from Z — > r + r with both rs decaying leptonically we 
reconstruct the t + t~ invariant mass, as explained below. 

5.2.1 Run I 



For the dilepton channel, our Run I analysis[28] proposed the following cuts: (1) 
$ T > 25 GeV;(2) at least one jet with p T > 15 GeV and \tj\ < 2; (3) an opposite sign 
pair with p T (£) > 8(5) GeV for £ = e(/i); and (4) an azimuthal dilepton angle 
satisfying 20° < /S.<p(£ + < 160°. We also defined a kinematical quantity called 
bigness, B (see eqn. (|5.8|))) and imposed the cut (5) B < 100 GeV. Our analysis, 
which ignored the important background Z — > t + t~, found a reach to 110 
GeV, < 90 GeV. We noted that in distinction to the 6-jet + lepton channel 

the signal here was not strongly attenuated in the limit — ► rn^, since a 6-tag 
was not required. 

This study was undertaken by the DO Collaboration on 75 ± 8 pb -1 of their 
Run I data, looking specifically in the dielectron channel ||, in which the signal is 
four times smaller than if both electrons and muons are counted. They used the 
RGSEARCH cut optimization method and imposed (1)# T > 22 GeV; (2) E T (j x ) > 
30 GeV; (3) p T (ei) > 16 GeV and p T (e 2 ) > 8 GeV; and (4) m(e+e - ) < 60 GeV. 
They also applied our bigness cut as (5) B < 90 GeV. Two events remained in 
the data set after these cuts, and the experimenters found too much Z — > t + t~ 
background to set a meaningful limit. Below, we outline a technique which may be 
used to control this background. 
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5.2.2 Run II 

The initial cuts for the dilepton channel are 
(i) # T > 25 GeV; 

(ii) at least one jet with p T > 15 GeV and \rj\ < 2; 

(m) an opposite-sign pair with p T (I) > 10 GeV and \q\ < 2; 

(w) 20° < A</>(£+,£-) < 160°; 

(u) m (£ + £~) not between 80 and 100 GeV for same-flavor 

After these cuts, less than 1 fb of Z — > = e, /x) remains, and we do not 

consider it any further. The other backgrounds, together with some representa- 



tive signal points, are displayed in Table 5.2. As shown in the table, the largest 
background after the Z veto cut (v) is Z — > t + t~ . Although the Z in such events 
cannot be reconstructed in as straightforward a manner as for Z — > £ + l~ , an indirect 
method is available. 

In the Z rest frame, the rs are highly relativistic, so that the lepton and 
neutrinos from r — » v T vt are strongly boosted along the r direction. We approximate 
the t + 3-momentum as r + ~ P+£ + where £ + is the observed £ + unit direction vector 
and P+ is the (unknown) magnitude of r + . Now, the Z transverse momentum 
all comes from its recoil against QCD radiation, so we can estimate it as Zt = 
— Shadrons £^r(had). From above, we can also write Z^ = (f + + r~) T = {P+l + + 
P_£~) T . Therefore after solving the two components of 

(p + l+ + p_r) x = -^^ T (had) 

for the two unknowns P+ and P_ we can reconstruct the r 4-momenta as 



r ± = W(P±^) 2 + m?, P±^). (5.7) 
79 





Figure 5.8: Distribution of m(r + ,r ) in the dilepton channel, after cut (v). The signal 
case is a 130 GeV t\ decaying to a 100 GeV W\. The vertical axes are in fb/bin. Only 
events for which to(t+,t~) > are shown; the percentage of such events is indicated in 
parentheses on the plots (so, 99.7% of the Z —> t + t~ events are plotted, but only 41% of 
the ti events). Dashed vertical lines indicate cut (vi). 
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Table 5.2: dilepton channel cross-sections, in fb, after cuts as discussed in the text. 
Signal levels required for detection are given in the "25%" (of background) and 
"5c" rows. The parenthesized 5a value in the last column is for an integrated 
luminosity of 25 fb -1 . For the signal cases, t\ and W\ masses are given in GeV, and 

m z x = m wJ 2 - 

Note that for this procedure to make sense for Z — > t + t~ both P + and P_ should 
be positive quantities, which is to say that neither r should be oppositely directed 
to its lepton. (This condition is also equivalent to requiring that Zt lie within the 
smaller angle of i\ and €^.) Distributions of m(r + ,T~) after cut (v) are displayed 
in Fig. |5.8| for those events with P+, P_ > 0. We make the cut 



{vi) m(r + ,T-) < 50 GeV or m(r + , t~) > 150 GeV 
where we define m(r + ,r _ ) = if either P± < 0. Essentially all of the Z — > t + t - 



events have P+,P_ > 0, while only about 40% of the signal events do. Table 5.2 
shows the results of applying cut (vi), where we see that only 2% of the Z — * t + t~ 
background survives, while about 5/6 of the signal is retained. 

Even after this cut, the signal is still below the 5<r level of observability. To 
achieve a further purification we note, as in the previous channel, that the presence 
in the background of two real Ws undergoing two-body decays W — > Iv will typically 
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B 



Figure 5.9: Distribution of B in the dilepton channel, after cut (vi). The signal case is a 
130 GeV t\ decaying to a 100 GeV Z\. The vertical axes are in fb/bin. Dashed lines at 
B = 120 indicate cut (vii). 
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result in harder lepton momenta. We also find, for the t\/W\/Z\ masses accessible in 
this channel, that the backgrounds generally have higher The background events 
are thus "bigger" than the signal events, and a convenient quantity summarizing 
these features is "bigness" 

B=p T (l + )+p T (r)+$ T (5.8) 



which was first introduced in our Run I work [28|. Figure 5.S shows distributions 
of B after cut (vi). We find that the discovery reach over the parameter space is 
maximized with the cut 

(vii) B < 120 GeV. 



Table 5.2 demonstrates the utility of this cut: only 20% of the ti and 50% of the 
WW backgrounds survive, while we keep 60-75% of the signal. 

The reach in the dilepton channel after cut (vii) is still substantially less than 
that of the 6-jet + lepton channel. The situation is improved a bit by considering 



the total jet transverse energy Jt = X) 



jets 



whose distributions are shown in 



Fig. 5.10. Jt for the W and Z backgrounds is rather lower than that of the signal, 



while Jp for tt is higher. We use the cut 

(viii) 25 < J T < 175 GeV 

which cuts 1/3 of the background at a cost of less than 1/6 in signal, as shown in 
Table K 



Figure 5.11 displays contours of the cross-sections after cut (viii) in the mi - 



m^ plane with m^ = m^j2 and B(W\ — > I) = B(W — > I). The total background 
is 49 fb. The 5<r discovery limit for 2 fb _1 is 25 fb, shown as a heavy solid line in 
the figure. The dashed line is at 12 fb = 25% of background. The reach in this 
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Figure 5.10: Distribution of Jt in the dilcpton channel, after cut (vii). The signal case is 
a 130 GeV i\ decaying to a 100 GeV Z\. The vertical axes are in fb/bin. Dashed lines show 
cut (vii) as 25 < J T < 175 GeV. 
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Figure 5.11: Cross-section contours, in fb, for the dilepton channel after cut (viii) as 
described in the text. The heavy solid line at 24.7 fb is the 5cr discovery limit for an 
integrated luminosity of 2 fb -1 . The dashed line is 12.1 fb of signal, which is 25% of 
background. The dot-dashed line shows the reach at 25 fb -1 after cut (ix) . 
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channel is significantly less than that of the fo-jet + lepton channel; t±s heavier than 
135 GeV are not accessible via this search. 

Relaxing the assumption B(W\ — ► t) « B(W — > I) 

Because the cross-section for dilepton events goes as the square of B(W\ — > £), 
substantial improvement is possible. The situation will continue to improve as the 



enhancement factor r of eqn. (5.6) discussed in the previous section increases to 
its limit r = 3, at which point the signal cross-section in this channel would go 
up by an order of magnitude. In such a case even just 2 fb _1 of data could find 
any stops with mi is m;, even if the chargino is very close to the kinematical 
limit rriyy < + mj. Also, should a stop be discovered in the 6-jet + lepton 
channel of the previous section, then a search in the dilepton channel might yield 
useful information about B(W% — ► £) which, together with basic facts about the W\ 
discernible from the discovery, could shed some light on other SUSY parameters. 

5.2.3 Run 11+ 

This channel is severely rate limited. We can reach the dashed 25% contour at 12.1 
fb in Fig. |5.11 by using cuts (i) — (viii) on an 8.3 fb _1 data sample, by eqn. 



This gives us stops out to is 150 GeV and charginos to i$ 120 GeV. We are 
background limited here and another cut is required to go further. For the 25 fb -1 
Run 11+ sample, we can eliminate much of the W and Z background by insisting 
on a tagged b quark: 

(ix) at least one tagged B. 

We did not impose a 6-tag requirement for the Run II cuts, since there is not enough 
signal at 2 fb -1 to sustain such a cut. Also, Fig. 5.11 shows that one pays a penalty 
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for the fo-tag, in that an unobservable region of small — is created where 



the b does not have enough energy to generate the tag. Table |5.2| exhibits the results 
of this cut. With (ix), the 25 fb" 1 5a discovery limit is extended to the dot-dashed 
line shown in Fig. |5.11| , which pushes the stop mass reach to just over 165 GeV, 
and increases the region where the SUSY signal can be identified in more than one 
channel. Detection in both channels helps to identify that the signal is indeed a 
stop, and can also yield useful information on other SUSY parameters relevant to 
branching fractions and masses. 
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Chapter 6 

Summary and Conclusions 



In this work we have studied two important stop decay modes, t\ — > cZ\ and 
t\ — > bW±, in the context of the Run II Tevatron experiments. We have identi- 
fied useful kinematical quantities for distinguishing the stop signal events from SM 
backgrounds. We found that these backgrounds could be controlled well enough to 
expose the stop over a significant range of the parameter space of many models, so 
that Tevatron experiments in the Main Injector era will probe models not accessible 
to LEP II. 



Our results are summarized in Fie. 6.1. Both the m? x m~ and mr x m~ 



planes are displayed, after the fashion of Fig. 2.6. As with that figure, the vertical 
axes for m-z and m~ will coincide in case m T -r, = 2m~ (indeed, the t\ — ► bW\ 
contours correspond to = 2m~^). 

For the t± — > cZ\ mode we show the discovery limits at 2 fb _1 and 25 fb _1 , 
for both the ^ T channel and the c-tagged channel. The irregular hatched region is 
excluded at the 95% confidence level by LEP II and Tevatron Run I experiments. 



Figure 6.1 shows the c-tagged channel channel to be clearly superior to the $ T 
channel for this decay mode. While the reach in t± mass is comparable for the two, 




Figure 6.1: Summary of 5a discovery limits for the light stop. The vertical axis on the 
left is neutralino mass, and that on the right is chargino mass. In case m~ = 2/71-^, these 

axes coincide as illustrated. Solid contours show the reach at 2 fb _1 ; dashed lines are for 25 
fb _1 . The irregularly shaped hatched region at the left is excluded at the 95% confidence 
level in the t\ — > cZ\ mode by LEP II and Tevatron Run I experiments. The hatched region 
under m~ = 93 is excluded by LEP II chargino mass limit. 
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c-tagging allows heavier m~ to be probed. This is chiefly due to the large jet -E+ 
requirement in the^ T + jets search which means that — m~ must be large to 
yield hard c-jets. In either channel, we may find stops with masses below 165 GeV in 
Run II. The c-tagged channel generally gives an extra 15 GeV of reach in Z% mass. 

For higher integrated luminosities, the stop mass may be pushed well past 
the top mass, and c-tagging with 25 fb _1 of data allows discovery of stops heavier 
than 200 GeV. In the Run 11+ scenario, c-tagging extends the Z\ mass reach by 20 
to 30 GeV over the ^ T channel. Note however that as the t\ becomes heavier the 
four-body decay discussed in Section |2.2| may have a substantial branching fraction. 

Results for the t\ — > bW\ decay mode are shown in the lower right portion 



of Fig. 6.1. The £ + 6-jet and dilepton channels are illustrated, under our usual 
assumptions rn^ = 2771-^ and B{W\ — > leptons) = B(W — ► leptons). The hatched 
region < 93 GeV has been excluded by LEP II at the 95% confidence level. 

For this mode, with VF-like chargino branching fractions, the dilepton signal 
is of scant use compared to the £ + 6-jet channel with just 2 fb _1 . However, when 
B(W\ — > leptons) S> B{W — ► leptons) the dilepton channel rapidly improves. A 
detection in either or both of these channels could, by differential analysis of the 
two, yield valuable information about other sectors of the MSSM, particularly if 
kinematical evidence could be used to constrain the SUSY masses. Even supposing 
B(W\ — > leptons) = B{W — > leptons), we could probe stops out past 165 GeV with 
the data of Run II, and push this another 20 GeV or so with the higher integrated 
luminosity of Run 11+ . 



Note that the Run 11+ discovery frontiers plotted in Fig. 6.1 for the i + b 
and dilepton channels coincide for < 165 GeV at the high end of the 

range. This reflects the fact that the £ + b signal is strongly attenuated in the limit 
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mi x — > as the available energy for a hard identifiable 6-jet decreases, while 

the dilepton suffers less in this limit. Since the statistical significance of the two 
channels is comparable in this region, the two could be statistically combined to 
gain a bit more reach in m~ . 

In conclusion, we have demonstrated that the Fermilab Tevatron can probe 
significant and interesting portions of the light stop parameter space with experi- 
ments conducted during the initial and extended Run II operating phases. Particu- 
larly important, the procedures we have developed here should allow the Tevatron to 
test most of the region ^ 165 GeV favored by electroweak baryogenesis within 
mSUGRA models with minimal low-energy particle content. 

The large top Yukawa coupling and the availability of substantial LR mixing 
due to the trilinear scalar coupling make the stop the lightest sfermion in many 
models, and it may, in principle, be the lightest charged supersymmetric particle. 
The fact that the Tevatron can subject this hypothesis to a strong test holds out 
the exciting possibility that supersymmetry may make its first appearance in the 
discovery of the stop. 



91 



Bibliography 



[1] Physics Analysis Workstation 2.12/10 Reference Manual, CERN Program Li- 
brary (1999) 

[2] S. Abachi et al. (DO Collaboration), "The DO Detector" Nucl. Instrum. Meth. 
A324, 185 (1995) 

[3] S. Abachi et al. (DO Collaboration), "Search for squarks and gluinos at yfs = 
1.8 TeV" Phys. Rev. Lett. 75, 618 (1995) 

[4] S. Abachi et al. (DO Collaboration), "Search for light top squarks in pp colli- 
sions at 1.8 TeV" Phys. Rev. Lett. 76, 2222 (1996) 

[5] S. Abachi et al. (DO Collaboration), "Direct measurement of the top quark 
mass" Phys. Rev. Lett. 79, 1197 (1997) 

[6] S. Abachi et al. (DO Collaboration), "Search for top squark pair production 
in the dielectron channel" Phys. Rev. D57, 589 (1998) 

[7] B. Abbot et al. (DO Collaboration), "Measurement of the top quark mass 
using dilepton events" Phys. Rev. Lett. 80, 2063 (1998) 



92 



[8] S. Abachi et al. (DO Collaboration), "Search for squarks and gluinos in events 
containing jets and a large imbalance in transverse momentum" Phys. Rev. 
Lett. 83, 4937 (1999) 

[9] B. Abbot et al. (DO Collaboration), "Direct measurement of top quark mass" 
Phys. Rev. D58, 052001 (1998) 

[10] G. Abbiendi et al. (OPAL Collaboration), "Search for scalar top and scalar 
bottom quarks at ^fs = 189 GeV at LEP" Phys. Lett. B456, 95 (1999) 

[11] G. Abbiendi et al. (OPAL Collaboration), "Search for chargino and neutralino 
production at ^fs = 189 GeV at LEP " CERN-EP/99-123, to appear in Eur. 
Phys. J. C (1999) 

[12] F. Abe et al. (CDF Collaboration), "The CDF detector: an overview" Nucl. 
Instrum. Meth. A271, 387 (1988) 

[13] F. Abe et al. (CDF Collaboration), "Measurement of the top quark mass" 
Phys. Rev. Lett. 79, 1992 (1997) 

[14] F. Abe et al. (CDF Collaboration), "Search for gluinos and squarks at the 
Fermilab Tevatron collider" Phys. Rev. D56, R1357 (1997) 

[15] F. Abe et al. (CDF Collaboration), "Measurement of the top quark mass with 
the Collider Detector at Fermilab" Phys. Rev. Lett. 80, 2779 (1998) 

[16] F. Abe et al. (CDF Collaboration), "Kinematics of it events at CDF" Phys. 
Rev. Lett. 80, 2767 (1998) 

[17] M. Acciarri et al. (L3 Collaboration), "Searches for scalar quarks in e + e~ 
interactions at ^/s = 189 GeV" Phys. Lett. B471, 308 (1999) 

93 



[18] M. Acciarri et al. (L3 Collaboration), "Search for charginos and neutralinos 
in e+e- collisions at y/s = 189 GeV" Phys. Lett. B472, 420 (2000) 

[19] T. Affolder et al. (CDF Collaboration), "Search for scalar top and scalar bot- 
tom quarks in pp collisions at yfs = 1.8 TeV" hep-ex/bf 9910049, to appear 
in Phys. Rev. Lett. (1999) 

[20] T. Affolder et al. (CDF Collaboration), "Search for scalar top quark produc- 
tion in pp collisions at = 1.8 TeV" FERMILAB-PUB-99/340 hep-ex/bf 
9912018 submitted to Phys. Rev. Lett. (1999) 

[21] S. Ahmed et al. (CLEO Collaboration), CLEO Report CONF 99-10 (1999) 

[22] U. Amaldi, W. de Boer and H. Furstenau, "Comparison of grand unified theo- 
ries with electroweak and strong coupling constants measured at LEP" Phys. 
Lett. B260, 447 (1991) 

[23] H. Baer, M. Diaz, P. Quintana and X. Tata, "Impact of physical principles at 
very high energy scales on the superparticle mass spectrum" JHEP 0004, 016 
(2000) 

[24] H. Baer, M. Drees, R. Godbole, J. Gunion and X. Tata, "Phenomenology of 
light stops at the Fermilab Tevatron" Phys. Rev. D44, 725 (1991) 

[25] H. Baer, K. Cheung and J. Gunion, "A heavy gluino as the lightest supersym- 
metric particle" Phys. Rev. D59, 075002 (1999) 

[26] H. Baer, F. Paige, S. Protopopescu and X. Tata, in Proceedings of the Work- 
shop on Physics at Current Accelerators and Supercolliders, edited by J. 
Hewett, A. White and D. Zeppenfeld, Argonne National Laboratory (1993) 

94 



[27] H. Baer, F. Paige, S. Protopopescu and X. Tata, "ISAJET 7.48: A monte 
carlo event generator for pp, pp, and e + e~ interactions" hep-ph/bf 0001086 
(2000) 

[28] H. Baer, J. Sender and X. Tata, "Search for top squarks at the Fermilab 
Tevatron collider" Phys. Rev. D50, 4517 (1994) 

[29] R. Barate et al. (ALEPH Collaboration) "Searches for sleptons and squarks 
in e + e~ collisions at 189 GeV" Phys. Lett. B469, 303 (1999) 

[30] R. Barate et al. (ALEPH Collaboration) "Search for charginos and neutralinos 
in e + e~ collisions at centre-of-mass energies near 183 GeV and constraints on 
the MSSM parameter space" Eur. Phys. J. Cll, 193 (1999) 

[31] R. Barbieri and G.F. Giudice, "6 — > 57 decay and supersymmetry" Phys. Lett. 
B309, 86 (1993) 

[32] V. Barger, M.S. Berger and P. Ohmann, "Implications of supersymmetric 
grand unification" MAD/PH/826 (1994) 

[33] V. Barger and R. Phillips, Collider Physics Addison Wesley Publishing (1988) 

[34] A. Bartl, H. Eberl, S. Kraml, W. Majerotto and W. Porod, "Scalar top quark 
production at p + pT colliders" hep-ph/bf 9712412 (1997) 

[35] A. Bartl, H. Eberl, S. Kraml, W. Majerotto and W. Porod, in Proceedings 
of the 2nd Joint ECFA/DESY Study on Physics and detectors for a linear 
electron-positron collider (2000) 

[36] W. Beenakker, M. Kramer, T. Plehn, M. Spira, P.M. Zerwas, "Stop production 
at hadron colliders" Nucl. Phys. B515, 3 (1998) 

95 



[37] E.L. Berger, B.W. Harris and Z. Sullivan, "Single top squark production via 
R parity violating super symmetric couplings in hadron collisions" Phys. Rev. 
Lett. 83, 4472 (1999) 

[38] C. Boehm, A. Djouadi and M. Drees, "Light scalar top quarks and supersym- 
metric dark matter" Phys. Rev. D61, 095006 (2000) 

[39] C. Boehm, A. Djouadi and Y. Mambrini, "Decays of the lightest top squark" 
Phys. Rev. D61, 095006 (2000) 

[40] A. Bouquet, J. Kaplan and C.A. Savoy, "Low-energy constraints on super- 
gravity parameters" Nucl. Phys. B262, 299 (1985) 

[41] CDF Collaboration, FERMILAB-Pub-96/390-E (1996) 

[42] F. de Campos et al., in Proceedings of the Workshop "Physics at Run II: 
Super symmetry /Higgs Fermilab, Feb. - Nov. 1998" (1998) 

[43] M. Carena, M. Quiros and C.E.M. Wagner, "Electroweak baryogenesis and 
Higgs and stop searches at LEP and the Tevatron" Nucl. Phys. B524, 3 
(1998) 

[44] L. Clavelli, "New possibilities for a light gluino" hep-ph/bf 9908342 (1999) 

[45] S. Coleman and J. Mandula, "All possible symmetries of the S matrix" Phys. 
Rev. 159, 1251 (1967) 

[46] DO Collaboration, FERMILAB Pub-96/357-E (1996) 

[47] A. Datta, M. Guchait and K.K. Jeong, "Stop squark search at Tevatron in 
the light slepton scenario" Int. J. Mod. Phys. A14, 2239 (1999) 



96 



[48] D. Delepine, J.-M. Gerard, R. Gonzalez Felipe and J. Weyers, "A light stop 
and electroweak baryogenesis" Phys. Lett. B386, 183 (1996) 

[49] S. Deser and B. Zumino, "Consistent supergravity" Phys. Lett. 62B, 335 
(1976) 

[50] D. Dicus and V. Mathur, Phys. Rev. D7, 3111 (1973) 

[51] R. Demina, J.D. Lykken, K.T. Matchev and A. Nomerotski, "Stop and sbot- 
tom searches in Run II of the Fermilab Tevatron" hep-ph/bf 9910275 (Sub- 
mitted to Phys. Rev. D) (1999) 

[52] S. Dimopoulos and D. Sutter, "The supersymmetric flavor problem" Nucl. 
Phys. B452, 496 (1995) 

[53] P.J. Dornan for the ALEPH Collaboration, "ALEPH status report" at LEPC 
March 7, 2000, Imperial College, London (2000) 

[54] C.B. Dover, T. Gaisser and G. Steigman, "Cosmological constraints on new 
stable hadrons" Phys. Rev. Lett. 42, 1117 (1979) 

[55] M. Drees and K. Hagiwara, "Supersymmetric contribution to the electroweak 
p parameter" Phys. Rev. D42, 1709 (1990) 

[56] M. Drees and S. Martin, in Electroweak Symmetry Breaking and New Physics 
at the TeV Scale edited by T. Barklow, S. Dawson, H. Haber and J. Siegrist, 
World Scientific (to be published) 

[57] F. Dyson, Sci. Am. sept. 64, 129 (1964) 

[58] J. Ellis and S. Rudaz, "Search for supersymmetry in toponium decays" Phys. 
Lett. 128B, 248 (1983) 

97 



[59] J. Erler and D.M. Pierce, "Bounds on super symmetry from electroweak pre- 
cision analysis" Nucl. Phys. B526, 53 (1998) 

[60] D. Freedman, S. Ferrara and P. van Nieuwenhuizen, "Progress toward a theory 
of supergravity" Phys. Rev. D13, 3214 (1976) 

[61] H. Georgi and S.L. Glashow, "Unity of all elementary particle forces" Phys. 
Rev. Lett. 32, 438 (1974) 

[62] M.T. Grisaru, W. Siegel and M. Rocek, "Improved methods for supergraphs" 
Nucl. Phys. B159, 429 (1979) 

[63] L. Girardello and M. T. Grisaru, "Soft Breaking Of Supersymmetry" Nucl. 
Phys. B194, 65 (1982) 

[64] R. Haag, J.T. Lopuszanski and M.F. Sohnius, "All possible generators of su- 
persymmetries of the S matrix" Nucl. Phys. B88, 257 (1975) 

[65] H. Haber and G. Kane, "The search for supersymmetry: probing physics 
beyond the standard model" Phys. Rept. 117, 75 (1985) 

[66] K. Hikasa and M. Kobayashi, "Light scalar top at e + e~ colliders" Phys. Rev. 
D36, 724 (1987) 

[67] M. Hosch, R.J. Oakes, K. Whisnant, Jin Min Yang, Bing-lin Young, X. Zhang, 
"Probing top quark decay into light stop in the super symmetric standard 
model at the upgraded Tevatron" Phys. Rev. D58, 034002 (1998) 

[68] L. Ibanez and C. Lopez, "N=l supergravity, the weak scale and the low-energy 
particle spectrum" Nucl. Phys. B233, 511 (1984) 



98 



[69] K. Inoue, A. Kakuto, H. Komatsu and H. Takeshita, "Aspects of grand unified 
models with softly broken super symmetry" Prog. Theor. Phys. 68, 927 (1982) 

[70] D.R.T. Jones, "Introduction to super symmetry" in S. Raby and T. Walker 
(eds.) The Building Blocks of Creation, Proceedings of the 1993 TASI, World 
Scientific (1994) 

[71] M. Jones, "Top mass measurements" talk given at 1998 Moriond Electroweak 
Conference (1998) 

[72] G. Jungman and M. Kamionkowski, "Supersymmetric dark matter" Phys. 
Rept. 267, 195 (1996) 

[73] S. Katsanevas and S. Melachroinos, SUSYGEN 

[74] V.A. Kuzmin, V.A. Rubakov and M.E. Shaposhnikov, "On the anomalous 
electroweak baryon number nonconservation in the early universe" Phys. Lett. 
B155, 36 (1985) 

[75] H.L. Lai et al, "Global QCD analysis and the CTEQ parton distributions" 
Phys. Rev. D51, 4763 (1995) 

[76] B. W. Lee, C. Quigg and H. B. Thacker, "Weak Interactions At Very High- 
Energies: The Role Of The Higgs Boson Mass" Phys. Rev. D16, 1519 (1977) 

[77] Chong Sheng Li, R.J. Oakes and Jin Min Yang, "One loop qcd corrections to 
top quark decay into neutralino and light stop" Phys. Rev. D54, 6883 (1996) 

[78] G. Marchesini, B.R. Webber, G. Abbiendi, I.G. Knowles, M.H. Seymour and 
L. Stanco, "HERWIG 6.1" Computer Phys. Commun. 67, 465 (1992) 



99 



[79] S. Mrenna, "SPYTHIA, a super symmetric extension of PYTHIA 5.7" Corn- 
put. Phys. Commun. 101 232 (1997) 

[80] G. Miinster, "The role of mathematics in contemporary theoretical physics" in 
E. Rudolf and I.-O. Stamatescu (eds.) Philosophy, Mathematics and Modern 
Physics, Springer Verlag (1994) 

[81] P. Nath, R. Arnowitt and A. Chamseddine, Applied N=l Supergravity, ICTP 
Series in Theoretical Physics Vol. 1, World Scientific (1984) 

[82] H.P. Nilles, "Supersymmetry, supergravity and particle physics" Phys. Rept. 
110, 1 (1984) 

[83] Y. Okada, "Light stop and the b -> sj process" Phys. Lett. B315, 119 (1993) 

[84] J.C. Pati and A. Salam, "Lepton number as the fourth color" Phys. Rev. Lett. 
32, 1083 (1974) 

[85] J.C. Pati, in Superstrings, Compositeness and Cosmology, edited by S. Gates 
and R. N. Mohapatra, World Scientific (1987) 

[86] W. Porod and T. Wdhrmann, "Higher order top squark decays" Phys. Rev. 
D55, 2907 (1997) 

[87] W. Porod, "More on higher order decays of the lighter top squark" Phys. Rev. 
D59, 5009 (1999) 

[88] A. Sakharov, "Violation of CP invariance, C asymmetry, and baryon asym- 
metry of the universe" JETP Lett. 5, 24 (1967) 

[89] J. Sender, "Top to stop branching fraction limits from the Fermilab Tevatron 
top experiments" Phys. Rev. D54, 3271 (1996) 

100 



[90] Lian You Shan, Shou Hua Zhu, "Top decays into lighter stop and gluino" 
hep-ph/bf 9811430 (1998) 

[91] A. Sopczak, "Higgs boson searches at LEP up to y/s = 202 GeV" presented at 
the Seventh International Symposium on Particles, Strings, and Cosmology, 
PASCOS-99, Granlibakken, USA, Dec 10-16, 1999 (2000) 

[92] L. Susskind, "The gauge hierarchy problem, technicolor, supersymmetry, and 
all that" Phys. Rept. 104, 181 (1984) 

[93] X. Tata, What is supersymmetry and how do we find it?, Lectures presented 
at the IX Jorge A. Swieca Summer School, Campos do Jordao, Brazil, Report 
no. UH-511-872-97 (1997) 

[94] Physics at the Tevatron Run II Plenary talk at Pheno 2000, Madison WI 
(2000) 

[95] S. Weinberg, The Quantum Theory of Fields, Volume I, Cambridge University 
Press (1995) 

[96] Eugene Wigner, Gruppentheorie und ihre Anwendung auf die Quanten- 
mechanik der Atomspektren, Braunschweig (1931) 

[97] Eugene Wigner, "The unreasonable effectiveness of mathematics in the natural 
sciences" in Symmetries and Reflections, Indiana University Press (1967) 



101 



